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Abstract 

The observable universe could be a 1+3-surface (the "brane") embedded in a 1+3+d- 
dimensional spacetime (the "bulk"), with Standard Model particles and fields trapped on the 
brane while gravity is free to access the bulk. At least one of the d extra spatial dimensions 
could be very large relative to the Planck scale, which lowers the fundamental gravity scale, 
possibly even down to the electroweak (~ TeV) level. This revolutionary picture arises in the 
framework of recent developments in M theory. The 1+10-dimensional M theory encompasses 
the known 1+9-dimensional superstring theories, and is widely considered to be a promising 
potential route to quantum gravity. At low energies, gravity is localized at the brane and 
general relativity is recovered, but at high energies gravity "leaks" into the bulk, behaving in 
a truly higher-dimensional way. This introduces significant changes to gravitational dynamics 
and perturbations, with interesting and potentially testable implications for high-energy astro- 
physics, black holes, and cosmology. Brane-world models offer a phenomenological way to test 
some of the novel predictions and corrections to general relativity that are implied by M theory. 
This review analyzes the geometry, dynamics and perturbations of simple brane-world mod- 
els for cosmology and astrophysics, mainly focusing on warped 5-dimensional brane-worlds 
based on the Randall-Sundrum models. We also cover the simplest brane-world models in 
which 4-dimensional gravity on the brane is modified at low energies - the 5-dimensional 
Dvali-Gabadadze-Porrati models. Then we discuss co-dimension two branes in 6-dimensional 
models. 
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Update (21 April 2010) 

Inserted subsections on numerical solutions in Sections [5] and [7] and five figures. Added new 
Subsection 12.21 on the Randall-Sundrum model in string theory. Added new Section [5] on Dvali- 
Gabadadze-Porrati models and nine figures. Added new Section [TO] on 6-dimensional models and 
one figure. More than 80 new references. 

Page 3: Extended a paragraph moved here from the end of Section [L2l 

Page II: Added Subsection 12.21 "RS model in string theory". 

Page 39: Added two paragraphs discussing the neglection of backrcaction due to metric per- 
turbations in the bulk. 

Page 56: Added Subsection 16.51 "Full numerical solutions". 

Page 66: Removed former Figure II: Damping of brane- world gravity waves on horizon re-entry 
due to massive mode generation. 

Page 67: Added Subsection 17.21 "Full numerical solutions". 

Page 76: Added Section H "DGP Models". 

Page 87: Added Section [TO] "6-Dimensional Models". 
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1 Introduction 



At high enough energies, Einstein's theory of general relativity breaks down, and will be superceded 
by a quantum gravity theory. The classical singularities predicted by general relativity in gravi- 
tational collapse and in the hot big bang will be removed by quantum gravity. But even below 
the fundamental energy scale that marks the transition to quantum gravity, significant corrections 
to general relativity will arise. These corrections could have a major impact on the behaviour 
of gravitational collapse, black holes, and the early universe, and they could leave a trace - a 
"smoking gun" - in various observations and experiments. Thus it is important to estimate these 
corrections and develop tests for detecting them or ruling them out. In this way, quantum gravity 
can begin to be subject to testing by astrophysical and cosmological observations. 

Developing a quantum theory of gravity and a unified theory of all the forces and particles of 
nature are the two main goals of current work in fundamental physics. There is as yet no generally 
accepted (pre-)quantum gravity theory. Two of the main contenders are M theory (for reviews 
see, e.g., |215l 13591 1380] ) and quantum geometry (loop quantum gravity; for reviews see, e.g., 
|368l 1412] ). It is important to explore the astrophysical and cosmological predictions of both these 
approaches. This review considers only models that arise within the framework of M theory. 

In this review, we focus on RS brane-worlds (mainly the RS 1-brane model) and their gen- 
eralizations, with the emphasis on geometry and gravitational dynamics (see [306 316j I27H 14271 
[HCTI 17701 ISMl mil SH HZ! for previous reviews with a broadly similar approach). Other 
reviews focus on string-theory aspects, e.g., jl481 13181 1971 1360] . or on particle physics aspects, e.g., 
[3571 1369] 1263] IT52] [75]. We also discuss the 5D DGP models, which modify general relativity 
at low energies, unlike the RS models; these models have become important examples in cos- 
mology for achieving late-time acceleration of the universe without dark energy. Finally, we give 
brief overviews of 6D models, in which the brane has co-dimension two, introducing very different 
features to the 5D case with co-dimension one branes. 



1.1 Heuristics of higher-dimensional gravity 

One of the fundamental aspects of string theory is the need for extra spatial dimensionf0. This 
revives the original higher-dimensional ideas of Kaluza and Klein in the 1920s, but in a new context 
of quantum gravity. An important consequence of extra dimensions is that the 

4-dimensional Planck scale Mp = M4 is no longer the fundamental scale, which is Mi+d, where 
d is the number of extra dimensions. This can be seen from the modification of the gravitational 
potential. For an Einstcin-Hilbert gravitational action we have 



^gravity - / d^x d^y ^^l^g V^+^^R - 2A4+ J , (1) 

(4+'^)G^B ^ ^'^'^Rab ~ \ ^'^'^R '^'^%AB = -A4+. + ^'^'^AB, (2) 

where = {x'^\y^, . . . and n^^d ^^^'^ gravitational coupling constant. 



f^l+d = 87rG4+d = (3) 

The static weak field limit of the field equations leads to the 4 + rf-dimensional Poisson equation, 
whose solution is the gravitational potential, 

- (4) 



^We do not consider timelike extra dimensions: see | 395| for an interesting example. 
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If the length scale of the extra dimensions is L, then on scales r < L, the potential is 4 + d- 
dimensional, V ^ 7--(i+''). By contrast, on scales large relative to L, where the extra dimensions 
do not contribute to variations in the potential, V behaves like a 4-dimensional potential, i.e., r L 
in the d extra dimensions, and V ^ L~''r~^. This means that the usual Planck scale becomes an 
effective coupling constant, describing gravity on scales much larger than the extra dimensions, 
and related to the fundamental scale via the volume of the extra dimensions: 

Ml^MlXiL''. (5) 

If the extra-dimensional volume is Planck scale, i.e., L ^ , then M^+d ^ ^Ip- But if the extra- 
dimensional volume is significantly above Planck scale, then the true fundamental scale M4+d can 
be much less than the effective scale Mp ^ 10^^ GeV. In this case, we understand the weakness of 
gravity as due to the fact that it "spreads" into extra dimensions and only a part of it is felt in 4 
dimensions. 

A lower limit on M^^d is given by null results in table-top experiments to test for deviations from 
Newton's law in 4 dimensions, V cx r^^ . These experiments currently [296] probe sub-millimetre 
scales, so that 

L < 10-1 mm - (10-1^ TcY)-^ Ah+d > lo(32-i5rf)/(<i+2) r^^y^ 

Stronger bounds for brane-worlds with compact fiat extra dimensions can be derived from null 
results in particle accelerators and in high-energy astrophysics |75) \87] \190\ I195j . 

1.2 Brane-worlds and M theory 

String theory thus incorporates the possibility that the fundamental scale is much less than the 
Planck scale felt in 4 dimensions. There are five distinct l4-9-dimensional superstring theories, 
all giving quantum theories of gravity. Discoveries in the mid-1990s of duality transformations 
that relate these superstring theories and the l-|-10-dimensional supergravity theory, led to the 
conjecture that all of these theories arise as different limits of a single theory, which has come to 
be known as M theory. The 11th dimension in M theory is related to the string coupling strength; 
the size of this dimension grows as the coupling becomes strong. At low energies, M theory can be 
approximated by l-|-10-dimensional supergravity. 

It was also discovered that p-branes, which are extended objects of higher dimension than strings 
(1-branes), play a fundamental role in the theory. In the weak coupling limit, p-branes {p > 1) 
become infinitely heavy, so that they do not appear in the perturbative theory. Of particular 
importance among p-branes are the D-branes, on which open strings can end. Roughly speaking, 
open strings, which describe the non-gravitational sector, are attached at their endpoints to branes, 
while the closed strings of the gravitational sector can move freely in the bulk. Classically, this is 
realised via the localization of matter and radiation fields on the brane, with gravity propagating 
in the bulk (see Figure [T]). 

In the Horava-Wittcn solution |204j . gauge fields of the standard model are confined on two 

1- |-9-branes located at the end points of an jZ^ orbifold, i.e., a circle folded on itself across a 
diameter. The 6 extra dimensions on the branes are compactified on a very small scale close to 
the fundamental scale, and their effect on the dynamics is felt through "moduli" fields, i.e., 5D 
scalar fields. A 5D realization of the Horava-Witten theory and the corresponding brane- world 
cosmology is given in |302[ 13031 1304) . 

These solutions can be thought of as effectively 5-dimensional, with an extra dimension that can 
be large relative to the fundamental scale. They provide the basis for the Randall-Sundrum (RS) 

2- brane models of 5-dimensional gravity |362| (see Figure [2]) . The single-brane Randall-Sundrum 
models [36 IJ with infinite extra dimension arise when the orbifold radius tends to infinity. The RS 
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Figure 1: Schematic of confinement of matter to the brane, while gravity propagates in the bulk 
(from [7S]). 
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models are not the only phenomenological realizations of M theory ideas. They were preceded by 
the Arkani-Hamed-Dimopoulos-Dvah (ADD) brane- world models [T51 [Til [T^ [SI 15701 li^nMlfTTI] . 

which put forward the idea that a large volume for the compact extra dimensions would lower the 
fundamental Planck scale, 

Mew - 1 TeV < A'/4+d < A/p ~ 10^6 TeV, (7) 

where Afp^ is the electroweak scale. If M4^ij is close to the lower limit in Equation (O, then this 
would address the long-standing "hierarchy" problem, i.e., why there is such a large gap between 
Afow and Afp. 

In the ADD models, more than one extra dimension is required for agreement with experiments, 
and there is "democracy" amongst the equivalent extra dimensions, which are typically flat. By 
contrast, the RS models have a "preferred" extra dimension, with other extra dimensions treated 
as ignorable (i.e., stabilized except at energies near the fundamental scale). Furthermore, this extra 
dimension is curved or "warped" rather than flat: The bulk is a portion of anti-de Sitter (AdSs) 
spacetime. As in the Horava-Wittcn solutions, the RS brancs arc Z2-symmetric (mirror symmetry), 
and have a tension, which serves to counter the influence of the negative bulk cosmological constant 
on the branc. This also means that the self-gravity of the brancs is incorporated in the RS models. 
The novel feature of the RS models compared to previous higher-dimensional models is that the 
observable 3 dimensions arc protected from the large extra dimension (at low energies) by curvature 
rather than straightforward compactification. 

The RS brane-worlds and their generalizations (to include matter on the brane, scalar fields 
in the bulk, etc.) provide phenomenological models that reflect at least some of the features of 
M theory, and that bring exciting new geometric and particle physics ideas into play. The RS 
models also provide a framework for exploring holographic ideas that have emerged in M theory. 
Roughly speaking, holography suggests that higher-dimensional gravitational dynamics may be 
determined from knowledge of the fields on a lower-dimensional boundary. The AdS/CFT corre- 
spondence is an example, in which the classical dynamics of the higher-dimensional gravitational 
field arc equivalent to the quantum dynamics of a conformal field theory (CFT) on the boundary. 
The RS model with its AdSs metric satisfies this correspondence to lowest perturbative order |130j 
(see also [MSl[3I8l[T94l[389||Ml[Mll350l2^ for the AdS/CFT correspondence in a cosmological 
context). 

Before turning to a more detailed analysis of RS brane-worlds. We discuss the notion of Kaluza- 
Klcin (KK) modes of the graviton. 



1.3 Heuristics of KK modes 

The dilution of gravity via extra dimensions not only weakens gravity on the branc, it also extends 
the range of graviton modes felt on the brane beyond the massless mode of 4-dimcnsional gravity. 
For simplicity, consider a flat brane with one flat extra dimension, compactificd through the iden- 
tification J/ -H- y + 27rnL, where n = 0, 1, 2, . . . . The perturbative 5D graviton amplitude can be 
Fourier expanded as 

/(x^y) = 5]e-«/V«(x'^), (8) 

n 

where are the amplitudes of the KK modes, i.e., the effective 4D modes of the 5D graviton. 
To see that these KK modes are massive from the brane viewpoint, we start from the 5D wave 
equation that the massless 5D field / satisfies (in a suitable gauge): 

(^ta/ = => Df + dlf = 0. (9) 
It follows that the KK modes satisfy a 4D Klein-Gordon equation with an effective 4D mass m„, 

□/« = ml /„, m„ = ^. (10) 
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The massless mode fo is the usual 4D graviton mode. But there is a tower of massive modes, 
L~^, 2L~^, . . . , which imprint the effect of the 5D gravitational field on the 4D brane. Compactness 
of the extra dimension leads to discreteness of the spectrum. For an infinite extra dimension, 
L — > oo, the separation between the modes disappears and the tower forms a continuous spectrum. 
In this case, the coupling of the KK modes to matter must be very weak in order to avoid exciting 
the lightest massive modes with m > 0. 

From a geometric viewpoint, the KK modes can also be understood via the fact that the 
projection of the null graviton 5-momentum ^^^pA onto the brane is timelike. If the unit normal to 
the brane is n^, then the induced metric on the brane is 

gAB = ''%AB - rLAUB, ^^^gABn-^n^ ^ I, QABn^ ^ 0, (H) 

and the 5-momcntum may be decomposed as 

^%A = mnA+PA, PAn'^ = 0, m^^^An"^, (12) 

where pA = gAB '^^p^ is the projection along the brane, depending on the orientation of the 
5-momentum relative to the brane. The effective 4- momentum of the 5D graviton is thus pA- 
Expanding ^^^gAB ^^'^P^ ^^V^ ~ 0, we find that 

- -m''. (13) 

It follows that the 5D graviton has an effective mass m on the brane. The usual 4D graviton 
corresponds to the zero mode, m = 0, when ^^^pA is tangent to the brane. 

The extra dimensions lead to new scalar and vector degrees of freedom on the brane. In 5D, 
the spin-2 graviton is represented by a metric perturbation ^^^/i^s that is transverse traceless: 

^^\^A = = as '^^a''- (14) 

In a suitable gauge, '^^^Hab contains a 3D transverse traceless perturbation /ly , a 3D transverse 
vector perturbation E^, and a scalar perturbation /3, which each satisfy the 5D wave equation [131| : 

(15) 
(16) 
(17) 

= 0. (18) 

The other components of ^^^Hab are determined via constraints once these wave equations are 
solved. The 5 degrees of freedom (polarizations) in the 5D graviton arc thus split into 2 {hij) + 2 
(Si) +1 (/?) degrees of freedom in 4D. On the brane, the 5D graviton field is felt as 

• a 4D spin-2 graviton hij (2 polarizations), 

• a 4D spin-1 gravi- vector (gravi-photon) Sj (2 polarizations), and 

• a 4D spin-0 gravi-scalar (3. 

The massive modes of the 5D graviton are represented via massive modes in all three of these fields 
on the brane. The standard 4D graviton corresponds to the massless zero-mode of hij. 

In the general case of d extra dimensions, the number of degrees of freedom in the graviton 
follows from the irreducible tensor representations of the isometry group as ^{d + l){d + A) . 
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2 Randall— Sundrum Brane- Worlds 



RS brane-worlds do not rely on compactification to localize gravity at the brane, but on the 
curvature of the bulk (sometimes called "warped compactification"). What prevents gravity from 
'leaking' into the extra dimension at low energies is a negative bulk cosmological constant, 

A5--| = -6/i^ (19) 

where i is the curvature radius of AdSs and ^ is the corresponding energy scale. The curvature 
radius determines the magnitude of the Ricmann tensor: 



(20) 



The bulk cosmological constant acts to "squeeze" the gravitational field closer to the brane. We 
can see this clearly in Gaussian normal coordinates = [x^^^y) based on the brane at y = 0, for 
which the AdSs metric takes the form 

^^W- = e-^^y^/'^Tj^^dx^dx" + dy^, (21) 

with r]fj,^ being the Minkowski metric. The exponential warp factor reflects the confining role of 
the bulk cosmological constant. The Z2-symmetry about the brane at y = is incorporated via 
the \y\ term. In the bulk, this metric is a solution of the 5D Einstein equations, 

(5)Gab = -As (^Ws, (22) 

i.e., ^^^Tab = in Equation ([2]). The brane is a fiat Minkowski spacetime, gABix"", 0) = rjf^^S^^ aS"^ b, 
with self-gravity in the form of brane tension. One can also use Poincare coordinates, which bring 
the metric into manifestly conformally flat form, 

(5)ds2 ^ _ ^rj^^dx'^dx'' + dz^] , (23) 

where z = le^l^ . 



The two RS models arc distinguished as follows: 

RS 2-brane: There are two branes in this model |362| . at y = and y = L, with Z'2-symmetry 
identiflcations 

y^-y, y + L^ L-y. (24) 
The branes have equal and opposite tensions ±A, where 

The positive-tension brane has fundamental scale M5 and is "hidden" . Standard model flelds 
are conflned on the negative tension (or "visible") brane. Because of the exponential warping 
factor, the effective Planck scale on the visible brane at y = L is given by 

So the RS 2-brane model gives a new approach to the hierarchy problem: even if M5 ~ ~ 
TeV, we can recover Mp ~ 10^^ TeV by choosing L/t large enough. Because of the finite 



-,2L/ 



(26) 
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separation between the branes, the KK spectrum is discrete. Furthermore, at low energies 
gravity on the branes becomes Brans-Dicke-Uke, with the sign of the Brans-Dicke parameter 
equal to the sign of the branc tension jl56) . In order to recover 4D general relativity at low 
energies, a mechanism is required to stabilize the inter-brane distance, which corresponds to 
a scalar field degree of freedom known as the radion |1751 141H I338| I285j . 

RS 1-brane: In this model |361) . there is only one, positive tension, brane. It may be thought of 
as arising from sending the negative tension brane off to infinity, L — )■ oo. Then the energy 
scales are related via 

„ Ml 

Ml = -^. (27) 

The infinite extra dimension makes a finite contribution to the 5D volume because of the 
warp factor: 

j d^X^^g ^2 J d^x J^dye-^y/^ = ^J d^x. (28) 

Thus the effective size of the extra dimension probed by the 5D graviton is ^. 

We will concentrate mainly on RS 1-branc from now on, referring to RS 2-branc occasionally. The 
RS 1-branc models arc in some sense the most simple and geometrically appealing form of a brane- 
world model, while at the same time providing a framework for the AdS / CFT correspondence |1301 
[3451 13781 [T941 [389] l393l [292] l350l [TST] . The RS 2-brane introduce the added complication of radion 
stabilization, as well as possible complications arising from negative tension. However, they remain 
important and will occasionally be discussed. 

2.1 KK modes in RS 1-brane 

In RS 1-brane, the negative A5 is offset by the positive brane tension A. The fine-tuning in 
Equation p5|) ensures that there is a zero effective cosmological constant on the brane, so that the 
brane has the induced geometry of Minkowski spacctime. To see how gravity is localized at low 
energies, we consider the 5D graviton perturbations of the metric |36H 11561 1172| [123] . 

'^%AB ^ (^Wb + e-^\y\l' ('^hAB, ^'^hAy = = (5)/^% = (^V^. (29) 

(see Figure [3]). This is the RS gauge, which is different from the gauge used in Equation ^5]) . but 
which also has no remaining gauge freedom. The 5 polarizations of the 5D graviton are contained 
in the 5 independent components of ^'^-'/i^jy in the RS gauge. 

We split the amplitude / of ^^^hj^s into 3D Fourier modes, and the linearized 5D Einstein 
equations lead to the wave equation [y > 0) 



2 f _ p-2y/£ 
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(30) 



Separability means that we can write 

/(i,2/) = ^^™(i)/™(y), (31) 

m 

and the wave equation reduces to 

+ {m^ + fc2)(^,„ = 0, (32) 
/™-7/m + e''^'™'/™-0. (33) 
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Figure 3: Gravitational field of a small point particle on the brane in RS gauge. (Figure taken 
from [156].) 



The zero mode solution is 



and the m > solutions are 

ipmit) ~ Am+ cxp + fc2 + A^- cxp ^— m? + k"^ tj 



(34) 
(35) 



(36) 
(37) 



The boundary condition for the perturbations arises from the junction conditions, Equa- 
tion discussed below, and leads to f'{t,0) ~ 0, since the transverse traceless part of the 
perturbed energy-momentum tensor on the brane vanishes. This implies 



Co = 0, 

The zero mode is normalizable, since 



C„ 



Yi(m^) ™- 



< oo. 



(38) 



(39) 



cx r 



The 



Its contribution to the gravitational potential V ~ '■'^''/ioo gives the 4D result, V 
contribution of the massive KK modes sums to a correction of the 4D potential. For r <C one 
obtains 

y{r) « (40) 
which simply reflects the fact that the potential becomes truly 5D on small scales. For r ^ £, 

GM 



V{r) 



2£2 
^+3^ 



(41) 
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which gives the small correction to 4D gravity at low energies from extra-dimensional effects. These 
effects serve to slightly strengthen the gravitational field, as expected. 

Table-top tests of Newton's laws currently find no deviations down to O(10~^ mm), so that 
f < 0.1 mm in Equation (|4ip. Then by Equations (|25|) and (|27p . this leads to lower limits on the 
brane tension and the fundamental scale of the RS 1-brane model: 

A > (1 TcV)'^, Ms > 10^ TeV. (42) 

These limits do not apply to the 2-brane case. 

For the 1-branc model, the boundary condition. Equation p8|) . admits a continuous spectrum 
m > of KK modes. In the 2-brane model, f'{t, L) = must hold in addition to Equation (|38p . 
This leads to conditions on to, so that the KK spectrum is discrete: 

m„ = ^ e-^/^ where Ji(x„) = Y,ix„). (43) 

The limit Equation (|^^ indicates that there are no observable collider, i.e., O(TeV), signatures 
for the RS 1-branc model. The 2-brane model by contrast, for suitable choice of L and £ so 
that nil = ©(TeV), does predict collider signatures that arc distinct from those of the ADD 
models [1901 [TM] . 



2.2 RS model in string theory 

There have been interesting developments in realizing the Randall-Sundrum brane- world model in 
string theory. A concrete example was found in type IIB supergravity |171) . The bosonic part of 
the 10-dimensional action is given by (see for example j358| ) 



2 



f?n 1 . , F, 



e*C(4) A G(3) A G(3) + 5ioc ■ (44) 



Here gs denotes the string metric. We have also defined the combined three-flux, G(3) — F(3) — 
tHi^2)i where r = C(o) + ie^'^, and 

^(5) = J^(5) - ^^2) A i7(3) + ^5(2) A F(3) . (45) 

The term Sioc is the action of localized objects, such as brancs. 

There is a 10-dimensional solution similar to the one in the RS model 

dslo = h{y)-^r]^^dx''dx'' + h(y)^^ndy"'dy" , (46) 

where .x^ are four-dimensional coordinates and arc coordinates on the 6-dimensional compact 
manifold Mq. If there are N coincident D3-brancs on the manifold, the 10-dimensional spacetime 
is described by AdS^ x where is a five-dimensional Einstein manifold. The warp factor is 
given by |421) 



Hr)-[^-j , (47) 
where r is the distance from the D3 branes in the (jmn metric and L is given by 

L = AT:agsNa'^. (48) 
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Here gs is the string coupling constant, a' is related to the string mass scale ms as a' = I/'ti^ and 
a depends on X^. 

Near the D3 brancs, the geometry is AdS^ x . At large r, the product structure of the 
10-dimensional metric breaks down and r ceases to be a good coordinate, and the 6D extra- 
dimensional space M.Q is not a product of a five sphere and one-dimensional spacetime. This gives 
a minimum distance tq. This cut-off acts as a TeV brane in the RS model. On the other hand, 
for r > Tniax, the AdS geometry smoothly glues into a Calabi-Yau compactification. Thus, the 
compact manifold plays the role of the Planck brane in the RS model, acting as a cut-off of the 
spacetime. The minimum warp factor /i(ro)~^/^ is determined by the flux associated with ^3 and 
H3 and thanks to the warping, it is possible to realize a large hierarchy jl71) 

(49) 

max ) 
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3 Covariant Approach to Brane- World Geometry and Dy- 
namics 



The RS models and the subsequent generahzation from a Minkowski brane to a Friedmann- 
Robertson-Walker (FRW) brane iM. .217, 221 jJl|333 [210] [1461 were derived as solutions 
in particular coordinates of the 5D Einstein equations, together with the junction conditions at 
the Z2-symmctric braneH. A broader perspective, with useful insights into the inter-play between 
4D and 5D effects, can be obtained via the covariant Shiromizu-Maeda-Sasaki approach |391j . in 
which the brane and bulk metrics remain general. The basic idea is to use the Gauss-Codazzi 
equations to project the 5D curvature along the brane. (The general formalism for relating the 
geometries of a spacetime and of hypersurfaces within that spacetime is given in [426] .) 
The 5D field equations determine the 5D curvature tensor; in the bulk, they are 

^'^Gab = - As ^%AB + 4 ^"^Tab, (50) 

where ^^^Tab represents any 5D energy- momentum of the gravitational sector (e.g., dilaton and 
moduli scalar fields, form fields). 

Let yhe a, Gaussian normal coordinate orthogonal to the brane (which is at ?/ = without loss 
of generality), so that UAdX^ = dy, with being the unit normal. The 5D metric in terms of 
the induced metric on {y = const.} surfaces is locally given by 

^%AB = gAB + riATiB. ^^^ds^ = g^,^{x°',y)dx''dx'' + dy^. (51) 

The extrinsic curvature of {y — const.} surfaces describes the embedding of these surfaces. It can 
be defined via the Lie derivative or via the covariant derivative: 

Kab = ^£n9AB = 9A^ ^''^^cnB, (52) 

so that 

K[AB] = = Kabtl^ , (53) 

where square brackets denote anti-symmetrization. The Gauss equation gives the 4D curvature 
tensor in terms of the projection of the 5D curvature, with extrinsic curvature corrections: 

Rabcd = '^"'^REFGHgA^'gB^gc'^gD" + 2Ka[cKd\b. (54) 

and the Codazzi equation determines the change of Kab along {y — const.} via 

VbK^'a - VaK = ^""^Rbc gA^'n^, (55) 

where K — K^a- 

Some other useful projections of the 5D curvature are: 

'•''^Refgh gA^'gB^gc'^n" = 2V[aKb]c, (56) 
^'^Refgh gA^'n^gB^'n" = -E^Kab + KacK"" b, (57) 
^"•^RcD gA^gB^" = Rab - £„Aab - KKab + 2KacK^b- (58) 

The 5D curvature tensor has Weyl (tracefree) and Ricci parts: 

^'^Rabcd ^ ^'^Cacbd + I ((-^Wic ^'^Rd]b - ^%B[C ^'^RoiA) - l^%A[c ^%D]B ^'^R- (59) 

^ Brane- worlds without Z2 symmetry have also been considered; see e.g., | 228| for recent work. 
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3.1 Field equations on the brane 

Using Equations ([5(1)) and ([M)) . it follows that 

1 . 2 " 



(5) 



rp A B 

J ab9h .91/ 

1 



where ^^^T = ^^^T^a, and where 

is the projection of the bulk Weyl tensor orthogonal to n^. This tensor satisfies 

Sabtt-^ = = £[AB] ~ £a^, 



(60) 



(61) 



(62) 



by virtue of the Weyl tensor symmetries. Evaluating Equation (|5(7)) on the brane (strictly, as 
y — ^ ±0, since £ab is not defined on the brane |391j ) will give the field equations on the brane. 

First, we need to determine K^^, at the brane from the junction conditions. The total energy- 
momentum tensor on the brane is 



T. 



(63) 



where T^^ is the energy-momentum tensor of particles and fields confined to the brane (so that 
Tabit-^ = 0). The 5D field equations, including explicitly the contribution of the brane, are then 



(5)g 



AB 



5 VAB 



4 



'Iab + J^ab 



Siy) 



(64) 



Here the delta function enforces in the classical theory the string theory idea that Standard Model 
fields are confined to the brane. This is not a gravitational confinement, since there is in general 
a nonzero acceleration of particles normal to the brane j305j . 

Integrating Equation (|64p along the extra dimension from y = —e to y = -l-e, and taking the 
limit e — )- 0, leads to the Israel-Darmois junction conditions at the brane, 



day - 3 



0, 



Kt,. - K„„ = -K 



Tnbranc _T-'branc 



(65) 

(66) 



where T 



bra 



gfivrpbianc ^ Thc Z2 Symmetry means that when you approach the brane from one 



side and go through it, you emerge into a bulk that looks the same, but with the normal reversed, 
n"^ — 7^ —n-^. Then Equation (|52|) implies that 



K 



(67) 



so that wc can use thc junction condition Equation ()66p to determine the extrinsic curvature on 
the brane: 



K„ 



1 



Tt,u + ^ (A - T) .g^^ 



(68) 



where T = T*^^, where we have dropped the (+), and where we evaluate quantities on the brane 
by taking the limit y — > +0. 
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Finally we arrive at the induced field equations on the brane, by substituting Equation (|68p 
into Equation 

^2 ^2 

G^y = -A.9p„ + n'^Tfj^ty + 6— - + 4— J'^t.. (69) 

The 4D gravitational constant is an effective coupling constant inherited from the fundamental 
coupling constant, and the 4D cosmological constant is nonzero when the RS balance between the 
bulk cosmological constant and the brane tension is broken: 

k^ = kI^ ^X4, (70) 

A = i [As + K^A] . (71) 

The first correction term relative to Einstein's theory is quadratic in the energy-momentum 
tensor, arising from the extrinsic curvature terms in the projected Einstein tensor: 

S^. = ^TT^,. - It^c.T% + ^g^. [ST^^T^P - T^] . (72) 

The second correction term is the projected Weyl term. The last correction term on the right of 
Equation ([69]) . which generalizes the field equations in |391) . is 



4 



g^Liy, (73) 



where ^^-'Tab describes any stresses in the bulk apart from the cosmological constant (see |313| for 
the case of a scalar field) . 

What about the conservation equations? Using Equations (f50|) . ([55|) and (|68p . one obtains 

V^T^, = -2 (^^T^snV^- (74) 

Thus in general there is exchange of energy-momentum between the bulk and the brane. From 
now on, we will assume that 

(5)Tab=0 => J-,,, =0, (75) 

so that 

^^^Gab = -As ^%AB in the bulk, (76) 

G^„ = -Ag^^ + K^T^^ + 6—5^,. - on the brane. (77) 

One then recovers from Equation (j74p the standard 4D conservation equations, 

VT^, = 0. (78) 

This means that there is no exchange of energy-momentum between the bulk and the brane; their 
interaction is purely gravitational. Then the 4D contracted Bianchi identities (VG^^ = 0), applied 
to Equation ([69]) . lead to 

V^f^, = ^ V^5^„ (79) 

which shows qualitatively how 1 + 3 spacetime variations in the matter-radiation on the brane can 
source KK modes. 

The induced field equations ([77)1 show two key modifications to the standard 4D Einstein field 
equations arising from extra-dimensional effects: 
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• S^i, ^ [T^uf' is the high-energy correction term, which is negligible for p ^ A, but dominant 
for p > A: 

^ ^ A- ^^^^ 

• is the projection of the bulk Weyl tensor on the brane, and encodes corrections from 5D 
graviton effects (the KK modes in the linearized case). From the brane-observer viewpoint, 
the energy-momentum corrections in S^^, are local, whereas the KK corrections in are 
nonlocal, since they incorporate 5D gravity wave modes. These nonlocal corrections cannot 
be determined purely from data on the brane. In the perturbative analysis of RS 1-brane 
which leads to the corrections in the gravitational potential. Equation (|^T|) . the KK modes 
that generate this correction are responsible for a nonzero Efi^] this term is what carries the 
modification to the weak-field field equations. The 9 independent components in the tracefree 
Efj^v are reduced to 5 degrees of freedom by Equation ((7^ ; these arise from the 5 polarizations 
of the 5D graviton. Note that the covariant formalism applies also to the two-brane case. In 
that case, the gravitational influence of the second brane is felt via its contribution to £fj,i/. 

3.2 5-dimensional equations and the initial- value problem 

The effective field equations are not a closed system. One needs to supplement them by 5D 
equations governing £^^, which are obtained from the 5D Einstein and Bianchi equations. This 
leads to the coupled system |377j 

£n-f^pi/ = K^c,K°'i, — Efj^i, — -As^p^, (81) 

S'nEtiu = V"SQ(py) + -A5 {Kf^„ — gfj.uK) + K^-^Rpiavp + ?>K°'(^^£y)a -~ K£^^ 

+ {K^aK^p - K^pK^,) K^^, (82) 

£n'S^j/Q ~ —'2^[^J.£v\a + Ka^B^yfi — 2Bap[^K^]'^ , (83) 
£n-R^i/Q/3 = —2Rf^^y[aK fj]'' — ^iBafiv + ^ fiBjiau, (84) 

where the "magnetic" part of the bulk Weyl tensor, counterpart to the "electric" part £^^^ is 

6^^" = 9f.^gu''gcP ^"'^CABCon''. (85) 
These equations arc to be solved subject to the boundary conditions at the brane, 

V^£^, = 4V^5^., (86) 

B^ua = 2V[^A'j,]„ = -KgVf^ \ T^]a - ■7:gy\aT] , (87) 



brane- 



where A = B denotes ^|brano = ^|br 

The above equations have been used to develop a covariant analysis of the weak field |377| . 
They can also be used to develop a Taylor expansion of the metric about the brane. In Gaussian 
normal coordinates. Equation (|51|) . we have — d/dy. Then we find 



g^Lu{x,y) = gf,^ix,0) - Kg 



\y\ 



-£^, + \4 (t^o.T% + ^(A - T)T^J\ + 1 (^4{X - T)2 - As ) g,. 



y=o+ 

(88) 



16 



In a non-covariant approach based on a specific form of tlie bulk metric in particular coordinates, 
the 5D Bianchi identities would be avoided and the equivalent problem would be one of solving 
the 5D field equations, subject to suitable 5D initial conditions and to the boundary conditions 
Equation (|68p on the metric. The advantage of the covariant splitting of the field equations 
and Bianchi identities along and normal to the brane is the clear insight that it gives into the 
interplay between the 4D and 5D gravitational fields. The disadvantage is that the splitting is 
not well suited to dynamical evolution of the equations. Evolution off the timclike brane in the 
spacelike normal direction does not in general constitute a well-defined initial value problem jl2j . 
One needs to specify initial data on a 4D spacelike (or null) surface, with boundary conditions 
at the brane(s) ensuring a consistent evolution |332| 1208] . Clearly the evolution of the observed 
universe is dependent upon initial conditions which are inaccessible to brane-bound observers; this 
is simply another aspect of the fact that the brane dynamics is not determined by 4D but by 
5D equations. The initial conditions on a 4D surface could arise from models for creation of the 
5D universe [1551 12591 [TO] l44l 1382] . from dynamical attractor behaviour |337| or from suitable 
conditions (such as no incoming gravitational radiation) at the past Cauchy horizon if the bulk is 
asymptotically AdS. 



3.3 The brane viewpoint: A 1 + 3-covariant analysis 

Following |305) , a systematic analysis can be developed from the viewpoint of a brane-bound ob- 
server. (See also j229l .) The effects of bulk gravity arc conveyed, from a brane observer viewpoint, 
via the local (5^,^) and nonlocal (f^i/) corrections to Einstein's equations. (In the more general 
case, bulk effects on the brane are also carried by J^^i,, which describes any 5D fields.) The f^i/ 
term cannot in general be determined from data on the brane, and the 5D equations above (or 
their equivalent) need to be solved in order to find f^i^. 

The general form of the brane energy-momentum tensor for any matter fields (scalar fields, 
perfect fluids, kinetic gases, dissipative fluids, etc.), including a combination of different fields, can 
be covariantly given in terms of a chosen 4-velocity u'^ as 

Tf,i, ^ pu^u^ + ph^„ + TT^i, + q^u^ + q„u^. (89) 

Here p and p are the energy density and isotropic pressure, respectively, and 

hf,„ = g^i, + u^Ui, = '■^^g^iy - n^Ui, + u^u^ (90) 

projects into the comoving rest space orthogonal to on the brane. The momentum density and 
anisotropic stress obey 

9m = QifJ.), TT^J^ = T^iniy), (91) 

where angled brackets denote the spatially projected, symmetric, and tracefree part: 

W^p. (92) 
In an incrtial frame at any point on the brane, we have 

- (1, 0), V = diag(0, 1, 1, 1), = (0, V,), W^o = = Y. = W.^-W,,, (93) 

where i,j — \,2, 3. 

The tensor 5^1/, which carries local bulk effects onto the brane, may then be irreducibly de- 
composed as 

S^, = ^ [2p2 - 37r„^^"^] u^u, + ^ [2p2 + App + tt^^tt"^ - Aq^q'^] h^, 

- Y2 + ^rtTT^i' - ^7r„(^7r^)" + ^q(^q^) + ^W(^u^) - ^9"7i-„(^u^). (94) 
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This simplifies for a perfect fluid or minimally-coupled scalar field to 

= -^p [pUf_,u^ + {p + 2p) hf_,^] . (95) 

The traccfrec S^^i, carries nonlocal bulk effects onto the brane, and contributes an effective 
"dark" radiative energy-momentum on the brane, with energy density p£, pressure /3, momen- 
tum density q^, and anisotropic stress tt^^: 

-■j^f^!^ = P£ (^Uf,u^ + ^hf^i'^ + qluu + qtu^ + TT^^. (96) 

We can think of this as a KK or Weyl "fluid". The brane "feels" the bulk gravitational field 
through this effective fluid. More specifically: 

• The KK (or Weyl) anisotropic stress tt^^ incorporates the scalar or spin-0 ("Coulomb"), the 
vector (transverse) or spin-1 (gravimagnetic), and the tensor (transverse traceless) or spin-2 
(gravitational wave) 4D modes of the spin-2 5D graviton. 

• The KK momentum density incorporates spin-0 and spin-1 modes, and defines a velocity 

of the Weyl fluid relative to u^^ via = P£V^. 

• The KK energy density p£, often called the "dark radiation", incorporates the spin-0 mode. 

In special cases, symmetry will impose simplifications on this tensor. For example, it must 
vanish for a conformally flat bulk, including AdSs, 

"•^Ws conformally flat = 0. (97) 

The RS models have a Minkowski brane in an AdS^ bulk. This bulk is also compatible with an FRW 
brane. However, the most general vacuum bulk with a Friedmann brane is Schwarzschild-anti-de 
Sitter spacetime |3391 l45j . Then it follows from the FRW symmetries that 

Schwarzschild AdSsbulk, FRW brane: 9^ = = tt^^, (98) 

where p£ ~ only if the mass of the black hole in the bulk is zero. The presence of the bulk black 
hole generates via Coulomb effects the dark radiation on the brane. 

For a static spherically symmetric brane (e.g., the exterior of a static star or black hole) jl05j . 

static spherical brane: q^ ~ 0. (99) 

This condition also holds for a Bianchi I brane |309j . In these cases, tt^^ is not determined by 
the symmetries, but by the 5D field equations. By contrast, the symmetries of a Godel brane fix 

The branc-world corrections can conveniently be consolidated into an effective total energy 



density, pressure, momentum density, and anisotropic stress: 

Ptot = P + ^ (2p' - 37r^.^^'^) + p£, (100) 

Ptot =P+j^ (2p' + 4pp + TT^^n^^ - 4q^Q^) + ^, (101) 

q';' = q, + ^{2pq^~37r^.qn + qf,, (102) 

ttJTJ' ^m"- + ^ [^(P + ^P)^t^'^ - 37i-a(^7r^)" + 3(j(^g^)] + nf^^. (103) 
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These general expressions simplify in the case of a perfect fluid (or minimally coupled scalar field, 
or isotropic one-particle distribution function), i.e., for = = tt^ 



to 



Ptot 

Ptot 
tot 

^tot 



2A 



p 



.s 



(104) 

(105) 

(106) 
(107) 



Note that nonlocal bulk effects can contribute to effective imperfect fluid terms even when the 
matter on the brane has perfect fluid form: There is in general an effective momentum density and 
anisotropic stress induced on the brane by massive KK modes of the 5D graviton. 

The effective total equation of state and sound speed follow from Equations ()104p and (jlOSp as 



Wtot 



Ptot _ w + (1 + 2w)p/2X + ps/ip 

Ptot 



^2 ^ Ptot 
^tot — • 

Ptot 



l + p/2A + pe/p 
P + P , '^Ps 



p + X 9{p + p){l + p/X) 



1 



4p£ 



3(p + p)(l + p/A) 



(108) 
(109) 



where w = p/ p and ~ p/ P- At very high energies, i.e., p 3> A, we can generally neglect p£ (e.g. 
in an inflating cosmology), and the effective equation of state and sound speed are stiffened: 



Wtot « 2w + 1, 



1. 



(110) 



This can have important consequences in the early universe and during gravitational collapse. For 
example, in a very high-energy radiation era, w = 1/3, the effective cosmological equation of state 
is ultra-stiff: wtot ~ 5/3. In late-stage gravitational collapse of pressureless matter, w = 0, the 
effective equation of state is stiff, wtot ~ 1, and the effective pressure is nonzero and dynamically 
important. 



3.4 Conservation equations 

Conservation of T^j/ gives the standard general relativity energy and momentum conservation 
equations, in the general, nonlinear case: 

p + e{p + p)+ V^q^ + 2Ay,, + a^''^^, = 0, (111) 

4 - 

<i{^) + 3 ©9,. + V^p +{p + p)Af, + V^TT^^ + A'^TT^^ + (T^.^q'' - e^.^^Lo'^q"' = 0. (112) 

In these equations, an overdot denotes u'^Vi/, & = V^u^ is the volume expansion rate of the 
worldlines, — — ^(p) is their 4-acceleration, cr^^ = V^^tit^^ is their shear rate, and 

cjp = — i curlwp = W(^) is their vorticity rate. 
On a Friedmann brane, we get 

= Lj,, = <j^, = 0, e = 3H, (113) 
where H = a/a \s the Hubble rate. The covariant spatial curl is given by 

curlV^^ = e^„/jV"F'', curW^, = e^p^^V'^W^ (114) 



19 



where Sp^ap is the projection orthogonal to of the 4D brane alternating tensor, and is the 
projected part of the brane covariant derivative, defined by 

V^i^"-...;3 - (V^i^"-...^)^„ = . . . hp'\/,F^-,„s. (115) 

In a local inertial frame at a point on the brane, with u'^ ~ S'^q, wc have: ^ Aq = loq = (To/j — 
eoQ/3 — curl Vb = curl Wq^ , and 

V^F"-.../3 = . ..df)''V,F^-...k (local inertial frame), (116) 

where A: = 1, 2, 3. 

The absence of bulk source terms in the conservation equations is a consequence of having A5 
as the only 5D source in the bulk. For example, if there is a bulk scalar field, then there is energy- 
momentum exchange between the brane and bulk (in addition to the gravitational interaction) |313[ 
[2TI [3241 [Mil [2761 [1451 148]. 

Equation ((79)) may be called the "nonlocal conservation equation" . Projecting along gives 
the nonlocal energy conservation equation, which is a propagation equation for p£. In the general, 
nonlinear case, this gives 



4 
3 

1 

2A 



-2g^VpP + 3g^V"^^. + 3^^"V^g, + Sa'^'na^TT,'' - Sa'^'qpq, . (117) 



Projecting into the comoving rest space gives the nonlocal momentum conservation equation, which 
is a propagation equation for q^ : 

1 



4A 



-4(p + p)Vpp + 6{p + p)V'^7r^^ + Gq^TTi^p^) + 271^"^ V^p 



-A{p + p)Qqp + 2qpA''q„ + GA^q^q^ + 4q^cr"'^7r„^ 



The 1 + 3-covariant decomposition shows two key features: 



(118) 



• Inhomogeneous and anisotropic effects from the 4D matter-radiation distribution on the 
brane are a source for the 5D Weyl tensor, which nonlocally "backreacts" on the brane via 
its projection S^^. 

• There are evolution equations for the dark radiative (nonlocal, Weyl) energy (ps) and mo- 
mentum (g^) densities (carrying scalar and vector modes from bulk gravitons), but there is 
no evolution equation for the dark radiative anisotropic stress (tt^^) (carrying tensor, as well 
as scalar and vector, modes), which arises in both evolution equations. 

In particular cases, the Weyl anisotropic stress n^^ may drop out of the nonlocal conservation 
equations, i.e., when we can neglect ct^'^tt^^,, V^tt^^, and A'^tt^^. This is the case when we consider 
linearized perturbations about an FRW background (which remove the first and last of these terms) 
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and further when we can neglect gradient terms on large scales (which removes the second term) . 
This case is discussed in Section [HI But in general, and especially in astrophysical contexts, the 
'^fif terms cannot be neglected. Even when wc can neglect these terms, arises in the field 
equations on the brane. 

All of the matter source terms on the right of these two equations, except for the first term on 
the right of Equation ()118|) . are imperfect fluid terms, and most of these terms are quadratic in 
the imperfect quantities and tt^^. For a single perfect fluid or scalar field, only the V^p term 
on the right of Equation (|118p survives, but in realistic cosmological and astrophysical models, 
further terms will survive. For example, terms linear in tt^^ will carry the photon quadrupolc in 
cosmology or the shear viscous stress in stellar models. If there are two fluids (even if both fluids 
are perfect), then there will be a relative velocity Vp, generating a momentum density = pw^, 
which will serve to source nonlocal effects. 

In general, the 4 independent equations in Equations (|117p and (jllSp constrain 4 of the 9 
independent components of on the brane. What is missing is an evolution equation for tt^^, 
which has up to 5 independent components. These 5 degrees of freedom correspond to the 5 
polarizations of the 5D graviton. Thus in general, the projection of the 5-dimensional field equations 
onto the brane does not lead to a closed system, as expected, since there are bulk degrees of freedom 
whose impact on the brane cannot be predicted by brane observers. The KK anisotropic stress 
TT^j, encodes the nonlocality. 

In special cases the missing equation does not matter. For example, if tt^^ = by symmetry, as 
in the case of an FRW brane, then the evolution of is determined by Equations (|117p and ()118|) . 
If the brane is stationary (with Killing vector parallel to u^), then evolution equations are not 
needed for £^1^, although in general tt^^, will still be undetermined. However, small perturbations 
of these special cases will immediately restore the problem of missing information. 

If the matter on the brane has a perfect-fluid or scalar-field energy-momentum tensor, the local 
conservation equations (jllip and (|112p reduce to 

p + e{p+p)^o, (119) 

V^.P+{p+p)A^=0, (120) 
while the nonlocal conservation equations (|117p and (jllSp reduce to 

P£ + ^ep£ + V^ql + 2A^^ql + a^^nl, = 0, (121) 
^1) + \H + \^,Pe + \peA, + + A^i^% + a^^ql - e^^^qi = -i^±^V^p. (122) 

Equation shows that 

• if £^l_u ~ and the brane energy-momentum tensor has perfect fluid form, then the density p 
must be homogeneous, V^p = 0; 

• the converse does not hold, i.e., homogeneous density does not in general imply vanishing 

A simple example of the latter point is the FRW case: Equation (|122p is trivially satisfied, 
while Equation (|121[) becomes 

PS + AHp£ = 0. (123) 
This equation has the dark radiation solution 

P£^P£o{^y. (124) 
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If £fj,i, = 0, then the field equations on the brane form a closed system. Thus for perfect fluid 
branes with homogeneous density and = 0, the brane field equations form a consistent closed 
system. However, this is unstable to perturbations, and there is also no guarantee that the resulting 
brane metric can be embedded in a regular bulk. 

It also follows as a corollary that inhomogeneous density requires nonzero S^i^: 

V^p^O ^ f,,. ^0. (125) 

For example, stellar solutions on the brane necessarily have S^i, ^ in the stellar interior if it is 
non-uniform. Perturbed FRW models on the brane also must have Ef^^, ^ 0. Thus a nonzero f^^, 
and in particular a nonzero tt^^, is inevitable in realistic astrophysical and cosmological models. 



3.5 Propagation and constraint equations on the brane 

The propagation equations for the local and nonlocal energy density and momentum density are 
supplemented by further 1 + 3-covariant propagation and constraint equations for the kinematic 
quantities O, A^, w^, cr^^, and for the free gravitational field on the brane. The kinematic quantities 
govern the relative motion of neighbouring fundamental world-lines. The free gravitational field 
on the brane is given by the brane Weyl tensor Cfj^^ap. This splits into the gravito-clcctric and 
gravito-magnetic fields on the brane: 

Ef^i, ^ Cf^aupu^'u^ = E(^f,^-^, H^j,^ = ]^e^,afiC°''^u-tU^ = H(^^), (126) 
where E^i, is not to be confused with E^^. The Ricci identity for 

V[^V^]Mq = ^i?ai/^/3U'^, (127) 

and the Bianchi identities 

V'^C^.c./j = V[,, + \r9.]^ , (128) 

produce the fundamental evolution and constraint equations governing the above covariant quan- 
tities. The field equations are incorporated via the algebraic replacement of the Ricci tensor 
by the effective total energy-momentum tensor, according to Equation (|69p . The brane equations 
are derived directly from the standard general relativity versions by simply replacing the energy- 
momentum tensor terms p, . . . by ptot, ■ • ■ • For a general fluid source, the equations are given 
in [3051 1229] . In the case of a single perfect fluid or minimally-coupled scalar field, the equations 
reduce to the following nonlinear equations: 

• Generalized Raychaudhuri equation (expansion propagation): 

1 2 2 

e + -e2+a^,a^--2c.^t.^-V^A^-A^A^ + y(p+3p)-A = -^(2p+3p)^ -K^f- (129) 



Vorticity propagation: 



• Shear propagation: 
2 



2 1 

■'(m) + 3®'^^' + 2 curl - cr^^w" = 0. (130) 
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Gravito-electric propagation (Maxwell- Weyl E-dot equation) 



E(^^^) + QE^^ - curl H^^ + y (p + P)(r^,u - 2 A"£„^(^i/^) ~ icTat^^E^) " +uj°'eaf}{i,E^) 



K 



4p£C7f,, + 37rf^,,) + Ott^, + 3V(^.(75) + OA^^gf^ + 3a"(^^5^„ + 3w"£„^(^7r5/ .(132) 



• Gravito-magnctic propagation (Maxwell- Weyl H-dot equation): 

H{f_iL') + QHf^i, + curlE'p^ - 'icTaiti.H^)"' + uj°'eafj{^E[v)^ + 2A°'£ap[fj,E^)^ 



K 

• Vorticity constraint: 

• Shear constraint: 



curlTT^^ - 3w(^gJ^ + (7a{f,e^f^q^ 



Va^^ - curla;^ - -V^O + 2e ^.^c^lu" A" = -K^q^ 

• Gravito-magnctic constraint: 

curlcr^^ + V(^a;^) - i^^i, + 2A(^w^) = 0. 

• Gravito-electric divergence (Maxwell- Weyl div-E equation): 



(133) 
(134) 

(135) 
(136) 



^T^mP + '— ( 2V^P£ - 2e<7^ - 3V'^7r^, + 3a/ - 9e/"w,gf 



(137) 



3 A ■ ' 6 

• Gravito-magnetic divergence (Maxwell- Weyl div-H equation): 

^"Hf,, -K^{p + p)io^ + e^.c^T'^fsE''^ - SE^.uj" = 

K.^{p + p)^oj^ + y {8p£0j^ -Scurlqf^ - 3£/"a/^f^ - iT^f..^") ■ (138) 



Gauss-Codazzi equations on the brane (with = 0): 



i?^ + r©^ - cr^,ycr'"' - 2k'^p-2A = n^^ + 2n^ pE 



(139) 
(140) 



where i?^^, is the Ricci tensor for 3-surfaces orthogonal to on the brane, and = /i^'^i?^^. 

The standard 4D general relativity results are regained when A^^ ^ and E^i^ ~ 0, which sets 
all right hand sides to zero in Equations (dH EM [Ml [HI 1133 [131 [133 [131 [133 [133 [133 
I140p . Together with Equations (|1191 1120[ 11211 1122p . these equations govern the dynamics of the 
matter and gravitational fields on the brane, incorporating both the local, high-energy (quadratic 
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Figure 4: The evolution of the dimcnsionless shear parameter Jlshcar = /QH^ on a Bianchi I 
brane, for a.V= ^m?(fp model. The early and late-time expansion of the universe is isotropie, but 
the shear dominates during an intermediate anisotropic stage. (Figure taken from j309) .) 



energy-momentum) and nonloeal, KK (projected 5D Weyl) effects from the bulk. High-energy 
terms are proportional to p/A, and are significant only when p > X. The KK terms contain pg, q^, 
and TT^j,, with the latter two quantities introducing imperfect fluid effects, even when the matter 
has perfect fluid form. 

Bulk effects give rise to important new driving and source terms in the propagation and con- 
straint equations. The vorticity propagation and constraint, and the gravito-magnetic constraint 
have no direct bulk effects, but all other equations do. High-energy and KK energy density terms 
are driving terms in the propagation of the expansion Q. The spatial gradients of these terms 
provide sources for the gravito-electric field -E^j/. The KK anisotropic stress is a driving term in 
the propagation of shear cr^^ and the gravito-electric/gravito- magnetic fields, E and Hfj,^ respec- 
tively, and the KK momentum density is a source for shear and the gravito-magnetic field. The 4D 
Maxwell-Weyl equations show in detail the contribution to the 4D gravito-eleetromagnetic field 
on the brane, i.e., (E^i,,H^^), from the 5D Weyl field in the bulk. 

An interesting example of how high-energy effects can modify general relativistic dynamics 
arises in the analysis of isotropization of Bianchi spacetimes. For a Binachi type I brane. Equa- 
tion (fT40|) becomes [309] 

4" +5' 

if we neglect the dark radiation, where a and H are the average scale factor and expansion rate, 
and E is the shear constant. In general relativity, the shear term dominates as a — 0, but in the 
brane-world, the high-energy term will dominate if w > 0, so that the matter-dominated early 
universe is isotropie [309, 68 , 57] mH [373 [H iS] . This is illustrated in Figure IH 

Note that this conclusion is sensitive to the assumption that « 0, which by Equation (|12ip 
implies the restriction 

0. (142) 
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Relaxing this assumption can lead to non-isotropizing solutions [5| I94| 166) . 

The system of propagation and constraint equations, i.e., Equations (|119[ 11201 11211 I122p 
and (Hm HSni inil ISa inSl HSl inSl ISl ISZl \IM HSni [Uni , is exact and nonhnear, applicable 
to both cosmological and astrophysical modelling, including strong-gravity effects. In general the 
system of equations is not closed: There is no evolution equation for the KK anisotropic stress 
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4 Gravitational Collapse and Black Holes on the Brane 



The physics of brane- world compact objects and gravitational collapse is complicated by a number 
of factors, especially the confinement of matter to the brane, while the gravitational field can access 
the extra dimension, and the nonlocal (from the brane viewpoint) gravitational interaction between 
the brane and the bulk. Extra-dimensional effects mean that the 4D matching conditions on the 
brane, i.e., continuity of the induced metric and extrinsic curvature across the 2-surface boundary, 
are much more complicated to implement |16H I120| I425| I160j . High-energy corrections increase 
the effective density and pressure of stellar and collapsing matter. In particular this means that 
the effective pressure does not in general vanish at the boundary 2-surfacc, changing the nature 
of the 4D matching conditions on the brane. The nonlocal KK effects further complicate the 
matching problem on the brane, since they in general contribute to the effective radial pressure 
at the boundary 2-surface. Gravitational collapse inevitably produces energies high enough, i.e., 
p ^ A, to make these corrections significant. 

We expect that extra-dimensional effects will be negligible outside the high-energy, short- 
range regime. The corrections to the weak- field potential. Equation (|4ip . are at the second 
post-Newtonian (2PN) level |165| I211| . However, modifications to Hawking radiation may bring 
significant corrections even for solar-sized black holes, as discussed below. 

A vacuum on the brane, outside a star or black hole, satisfies the brane field equations 

= R\ = Q^E\, V'f^. = 0. (143) 

The Weyl term f^i, will carry an imprint of high-energy effects that source KK modes (as discussed 
above). This means that high-energy stars and the process of gravitational collapse will in general 
lead to deviations from the 4D general relativity problem. The weak-field limit for a static spherical 
source. Equation (j4T|) . shows that must be nonzero, since this is the term responsible for the 
corrections to the Newtonian potential. 



4.1 The black string 

The projected Weyl term vanishes in the simplest candidate for a black hole solution. This is 
obtained by assuming the exact Schwarzschild form for the induced brane metric and "stacking" 
it into the extra dimension [76], 

(5)ds2 ^ e-^\y\/^g^,dx^'dx'' + dy^, (144) 
~9,. = e^\y\/'g,, = -(1 - 2GM/r)de + ^ ^^'^^^^ + rHn\ (145) 

(Note that Equation (|144p is in fact a solution of the 5D field equations (|22p if g^i, is any 4D Einstein 
vacuum solution, i.e., if = 0, and this can be generalized to the case i?^^ = —Ag^^, [111 I20| .) 

Each {y = const.} surface is a 4D Schwarzschild spacetime, and there is a line singularity along 
r ~ for all y. This solution is known as the Schwarzschild black string, which is clearly not 
localized on the brane y = 0. Although ^^^Cabcd 0, the projection of the bulk Weyl tensor 
along the brane is zero, since there is no correction to the 4D gravitational potential: 

, , GM , , 

y{r)^— £^, = 0. (146) 

The violation of the perturbative corrections to the potential signals some kind of non-AdSs pathol- 
ogy in the bulk. Indeed, the 5D curvature is unbounded at the Cauchy horizon, as y — > oo [76) : 

^'^Rabcd (^)i?^^^^ = S + e^l^l/^. (147) 
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Furthermore, the black string is unstable to large-scale perturbations jl83) . 

Thus the "obvious" approach to finding a brane black hole fails. An alternative approach is 
to seek solutions of the branc field equations with nonzero E^i, |105| . Brane solutions of static 
black hole exteriors with 5D corrections to the Schwarzschild metric have been found |1051 11611 
I120| 1425) I225| [73l 1224) . but the bulk metric for these solutions has not been found. Numerical 
integration into the bulk, starting from static black hole solutions on the brane, is plagued with 
difiicuhies [3921178]. 



4.2 Taylor expansion into the bulk 

One can use a Taylor expansion, as in Equation (j88p , in order to probe properties of a static black 
hole on the brane |106) . (An alternative expansion scheme is discussed in [73].) For a vacuum 
brane metric. 



9ij,u{x, y) = gf,„{x, 0) - £^,u{x, 0+)y'^ - ^^^^.(a;, 0+)|y|' 



1 

12 



32 p 



y 

y=o+ 



(148) 



This shows in particular that the propagating effect of 5D gravity arises only at the fourth order 
of the expansion. For a static spherical metric on the brane. 



-F(r)df 



H{r) 



the projected Weyl term on the brane is given by 



F 



H' - 



1 - H 



rH 



F' l-H 
"f r 



pi jji 

Sge = -1 + H+-H[— + — 



(149) 

(150) 
(151) 
(152) 



These components allow one to evaluate the metric coefficients in Equation (|148p . For example, the 
area of the 5D horizon is determined by ggg] defining il}{r) as the deviation from a Schwarzschild 
form for H, i.e., 

2?71 

H{r)^l +V(r), (153) 



where m is constant, we find 



5ee(r,2/)=r2-7/.' I + -I2/IU2 + 



^' + -(1 + V'')(^V''-^/')' 



(154) 



This shows how ij) and its r-derivatives determine the change in area of the horizon along the extra 



dimension. For the black string '0 = 0, and we have ggg{r,y) 



horizon scale ^ i, we have from Equation (|4ip that 



3r3 



For a large black hole, with 



(155) 



This implies that ggg is decreasing as we move off the brane, consistent with a pancake- like shape of 
the horizon. However, note that the horizon shape is tubular in Gaussian normal coordinates |164l. 
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4.3 The "tidal charge" black hole 



The equations (|143p form a system of constraints on the brane in the stationary case, inchiding 
the static spherical case, for which 



e = = 

The nonlocal conservation equations E^i, 



1. 



reduce to 



where, by symmetry, 



(156) 



(157) 



(158) 



for some n£(r), with being the unit radial vector. The solution of the branc field equations 
requires the input of from the 5D solution. In the absence of a 5D solution, one can make an 
assumption about or p^i. to close the 4D equations. 

If we assume a metric on the brane of Schwarzschild-like form, i.e., H = F in Equation (|149p . 
then the general solution of the brane field equations is [105] 



f 



F = 1- 

r 

2G£Q 



2GM 2GiQ 



[u^Uiy - 2r^r^ + h 



(159) 
(160) 



where Q is a constant. It follows that the KK energy density and anisotropic stress scalar (defined 
via Equation (|158p ) arc given by 

Pe = ^.= (161) 

The solution (|159p has the form of the general relativity Reissner-Nordstrom solution, but 
there is no electric field on the brane. Instead, the nonlocal Coulomb effects imprinted by the bulk 
Weyl tensor have induced a "tidal" charge parameter Q, where Q = Q{M), since M is the source 
of the bulk Weyl field. We can think of the gravitational field of M being "reflected back" on the 
brane by the negative bulk cosmological constant |104| . If we impose the small-scale perturbative 
limit (r ^ £) in Equation (|40|) , we find that 



= -2M. 



(162) 



Negative Q is in accord with the intuitive idea that the tidal charge strengthens the gravitational 
field, since it arises from the source mass M on the branc. By contrast, in the Rcissncr-Nordstrom 
solution of general relativity, Q cx +(7^, where q is the electric charge, and this weakens the 
gravitational field. Negative tidal charge also preserves the spacelike nature of the singularity, and 
it means that there is only one horizon on the brane, outside the Schwarzschild horizon: 



rh = GM 



1 - 



2eQ 
GAP 



= GM 



GM 



(163) 



The tidal-charge black hole metric does not satisfy the far-field r^^ correction to the gravi- 
tational potential, as in Equation (|4ip . and therefore cannot describe the end-state of collapse. 
However, Equation p59p shows the correct 5D behaviour of the potential (oc r^^) at short dis- 
tances, so that the tidal-charge metric could be a good approximation in the strong- field regime 
for small black holes. 



28 



4.4 Realistic black holes 



Thus a simple brane-based approach, while giving useful insights, does not lead to a realistic black 
hole solution. There is no known solution representing a realistic black hole localized on the brane, 
which is stable and without naked singularity. This remains a key open question of nonlinear 
brane-world gravity. (Note that an exact solution is known for a black hole on a 1 + 2-brane in a 
4D bulk [142j . but this is a very special case.) Given the nonlocal nature of f^^. it is possible that 
the process of gravitational collapse itself leaves a signature in the black hole cnd-statc, in contrast 
with general relativity and its no-hair theorems. There arc contradictory indications about the 
nature of the realistic black hole solution on the brane: 

• Numerical simulations of highly relativistic static stars on the brane [430] indicate that 
general relativity remains a good approximation. 

• Exact analysis of Oppenheimer-Snyder collapse on the brane shows that the exterior is non- 
static [160] . and this is extended to general collapse by arguments based on a generalized 
AdS/CFT correspondence pSllTin] . 

The first result suggests that static black holes could exist as limits of increasingly compact static 
stars, but the second result and conjecture suggest otherwise. This remains an open question. 
More recent numerical evidence is also not conclusive, and it introduces further possible subtleties 
to do with the size of the black hole |264| 1433) . 

On very small scales relative to the AdSs curvature scale, r ^ £, the gravitational potential 
becomes 5D, as shown in Equation (|40|) . 



, G£M G5M 
^(0 ~ ^ - (164) 

In this regime, the black hole is so small that it does not "see" the brane, so that it is approximately 
a 5D Schwarzschild (static) solution. However, this is always an approximation because of the self- 
gravity of the brane (the situation is different in ADD-type brane-worlds where there is no brane 
tension). As the black hole size increases, the approximation breaks down. Nevertheless, one might 
expect that static solutions exist on sufficiently small scales. Numerical investigations appear to 
confirm this |264l 1433) : Static metrics satisfying the asymptotic AdSs boundary conditions are 
found if the horizon is small compared to ^, but no numerical convergence can be achieved close 
to I. The numerical instability that sets in may mask the fact that even the very small black holes 
are not strictly static. Or it may be that there is a transition from static to non-static behaviour. 
Or it may be that static black holes do exist on all scales. 

The 4D Schwarzschild metric cannot describe the final state of collapse, since it cannot incor- 
porate the 5D behaviour of the gravitational potential in the strong-field regime (the metric is 
incompatible with massive KK modes). A non-perturbative exterior solution should have nonzero 
£^1, in order to be compatible with massive KK modes in the strong-field regime. In the end-state 
of collapse, we expect an which goes to zero at large distances, recovering the Schwarzschild 
weak-field limit, but which grows at short range. Furthermore, £^,1, may carry a Weyl "fossil 
record" of the collapse process. 



4.5 Oppenheimer— Snyder collapse gives a non-static black hole 

The simplest scenario in which to analyze gravitational collapse is the Oppenheimer-Snyder model, 
i.e., collapsing homogeneous and isotropic dust |160) . The collapse region on the brane has an FRW 
metric, while the exterior vacuum has an unknown metric. In 4D general relativity, the exterior is 
a Schwarzschild spacetime; the dynamics of collapse leaves no imprint on the exterior. 
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The collapse region has the metric 

(1 + A:r2/4)2 ^ ' 

where the scale factor satisfies the modified Friedmann equation (see below), 

d^ 8ttG f p ps' 



The dust matter and the dark radiation evolve as 

P = Po(^) , P^=P£o(^) , (167) 

where oq is the epoch when the cloud started to collapse. The proper radius from the centre of 
the cloud is R{t) — 7'a(r)/(l + |fcr^). The collapsing boundary surface E is given in the interior 
comoving coordinates as a free-fall surface, i.e. r = tq = const., so that R^{t) = roa(T)/(l + jkr^). 
We can rewrite the modified Friedmann equation on the interior side of S as 

• , 2GM 3GAP Q , , 

where the "physical mass" M (total energy per proper star volume), the total "tidal charge" Q, 
and the "energy" per unit mass E are given by 

M = (169) 



3(l + ifcr2)3 

Now we assume that the exterior is static, and satisfies the standard 4D junction conditions. 
Then wc check whether this exterior is physical by imposing the modified Einstein equations p43p . 
We will find a contradiction. 

The standard 4D Darmois-Israel matching conditions, which we assume hold on the brane, 
require that the metric and the extrinsic curvature of S be continuous (there are no intrinsic 
stresses on E). The extrinsic curvature is continuous if the metric is continuous and if R is 
continuous. We therefore need to match the metrics and R across S. 

The most general static spherical metric that could match the interior metric on S is 



ds^ = -F{Rf 



^ 2Gm{R) 



R 



We need two conditions to determine the functions F{R) and m{R). Now E is a freely falling 
surface in both metrics, and the radial geodesic equation for the exterior metric gives R^ = —1 + 
2Gm{R) / R + E/F{R)'^, where is a constant and the dot denotes a proper time derivative, as 
above. Comparing this with Equation ()168|) gives one condition. The second condition is easier to 
derive if we change to null coordinates. The exterior static metric, with 

dR 

dv = dt 



(1 - 2Gm/R)F' 
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ds^ =^ -F'^ ( 1 — ) dw^ + 2F dv dR + dVL\ (173) 



R 

The interior Robertson- Walker metric takes the form |160| 



•2 [l~{k + a2)i?2/a2j ^^2 2T,,dvdR 

(1 - fci?2/a2) ^ ^1 - fci?2/a2 



where 



4'-' J rd - 1 + ifcr2' 

Comparing Equations (|173p and p74p on S gives the second condition. The two conditions together 
imply that F is a constant, which we can take as F{R) = 1 without loss of generality (choosing 
E = E + l), and that 

In the hmit A"^ = = Q, we recover the 4D Schwarzschild solution. In the general brane-world 
case, Equations (|172p and (|175p imply that the brane Ricci scalar is 



However, this contradicts the field equations ()143p . which require 

= 0- (177) 

It follows that a static exterior is only possible if M/X = 0, which is the 4D general relativity limit. 
In the brane-world, collapsing homogeneous and isotropic dust leads to a non-static exterior. Note 
that this no-go result does not require any assumptions on the nature of the bulk spacetime, which 
remains to be determined. 

Although the exterior metric is not determined (see [179] for a toy model), we know that its 
non-static nature arises from 

• 5D bulk graviton stresses, which transmit effects nonlocally from the interior to the exterior, 
and 

• the non- vanishing of the effective pressure at the boundary, which means that dynamical 
information from the interior can be conveyed outside via the 4D matching conditions. 

The result suggests that gravitational collapse on the brane may leave a signature in the exterior, 
dependent upon the dynamics of collapse, so that astrophysical black holes on the brane may in 
principle have KK "hair". It is possible that the non-static exterior will be transient, and will tend 
to a static geometry at late times, close to Schwarzschild at large distances. 



4.6 AdS/CFT and black holes on 1-brane RS-type models 

Oppenheimer-Snyder collapse is very special; in particular, it is homogeneous. One could argue 
that the non-static exterior arises because of the special nature of this model. However, the un- 
derlying reasons for non-static behaviour are not special to this model; on the contrary, the role 
of high-energy corrections and KK stresses will if anything be enhanced in a general, inhomoge- 
neous collapse. There is in fact independent heuristic support for this possibility, arising from the 
AdS/CFT correspondence. 
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The basic idea of the correspondence is that the classical dynamics of the AdSs gravitational 
field correspond to the quantum dynamics of a 4D conformal field theory on the brane. This 
correspondence holds at linear perturbative order [130) . so that the RS 1-brane infinite AdSs 
brane-world (without matter fields on the brane) is equivalently described by 4D general relativity 
coupled to conformal fields, 

G^, = ^nGTlf\ (178) 

According to a conjecture }409) . the correspondence holds also in the case where there is strong 
gravity on the brane, so that the classical dynamics of the bulk gravitational field of the brane 
black hole are equivalent to the dynamics of a quantum-corrected 4D black hole (in the dual 
CFT-plus-gravity description). In other words [409| 1140] : 



Quantum backreaction due to Hawking radiation in the 4D picture is described as classical 
dynamics in the 5D picture. 



• The black hole evaporates as a classical process in the 5D picture, and there is thus no 
stationary black hole solution in RS 1-brane. 

A further remarkable consequence of this conjecture is that Hawking evaporation is dramatically 
enhanced, due to the very large number of CFT modes of order (^/€p)^. The energy loss rate due 
to evaporation is 

^ nIt^], (179) 



where N is the number of light degrees of freedom. Using N ^ ^^/G, this gives an evaporation 
timescale [409] 

M \ ^ f 1 mm \ ^ 



Wap-^^j ^-^) yr. (180) 

A more detailed analysis |141j shows that this expression should be multiplied by a factor « 
100. Then the existence of stellar-mass black holes on long time scales places limits on the AdSs 
curvature scale that are more stringent than the table-top limit. Equation ([6]). The existence of 
black hole X-ray binaries implies 

£ < 10"^ mm, (181) 

already an order of magnitude improvement on the table-top limit. 

One can also relate the Oppenheimer-Snyder result to these considerations. In the AdS/CFT 
picture, the non- vanishing of the Ricci scalar. Equation (jl76p . arises from the trace of the Hawking 
CFT energy- momentum tensor, as in Equation (|178p . If we evaluate the Ricci scalar at the black 
hole horizon, R ~ 2GM, using A = 6Ml/M^, we find 

The CFT trace on the other hand is given by T^'^f') - NT^JM^, so that 

8^GT(^f') ^ (183) 

Thus the Oppenheimer-Snyder result is qualitatively consistent with the AdS/CFT picture. 

Clearly the black hole solution, and the collapse process that leads to it, have a far richer 
structure in the brane-world than in general relativity, and deserve further attention. In particular, 
two further topics are of interest: 
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• Primordial black holes in 1-brane RS-type cosmology have been investigated in \211\ 11861 
11851 13151 inU 1387] . High-energy effects in the early universe (see the next Section [S]) can 
significantly modify the evaporation and accretion processes, leading to a prolonged survival 
of these black holes. Such black holes evade the enhanced Hawking evaporation described 
above when they are formed, because they are much smaller than £. 

• Black holes will also be produced in particle collisions at energies > M5, possibly well below 
the Planck scale. In ADD brane- worlds, where Mi+d = ©(TeV) is not ruled out by current 
observations if d > 1, this raises the exciting prospect of observing black hole production 
signatures in the next-generation colliders and cosmic ray detectors (see [751 11701 1139) ). 
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5 Brane- World Cosmology: Dynamics 



A 1 + 4-dimensional spacetime with spatial 4-isotropy (4D splierical/plane/hyperbolic symmetry) 
has a natural foliation into the symmetry group orbits, which are 1 + 3-dimensional surfaces with 
3-isotropy and 3-homogeneity, i.e., FRW surfaces. In particular, the AdSs bulk of the RS brane- 
world, which admits a foliation into Minkowski surfaces, also admits an FRW foliation since it 
is 4-isotropic. Indeed this feature of 1-brane RS-type cosmological brane-worlds underlies the 
importance of the AdS/CFT correspondence in brane- world cosmology |345| I378| I194| I389| 13931 

The generalization of AdS5 that preserves 4-isotropy and solves the vacuum 5D Einstein equa- 
tion (|22p is Schwarzschild-AdSs, and this bulk therefore admits an FRW foliation. It follows 
that an FRW brane-world, the cosmological generalization of the RS brane-world, is a part of 
Schwarzschild-AdSs , with the 2'2-symmetric FRW brane at the boundary. (Note that FRW 
branes can also be embedded in non- vacuum generalizations, e.g., in Reissner-Nordstrom-AdSs 
and Vaidya-AdSs . ) 

In natural static coordinates, the bulk metric is 



(5)^52 ^ ^F{R)dT'^ 



F(R) 



l-Kr^ 



r^dn^ 



F{R) = K 



ie_ 



m 
R^' 



(184) 
(185) 



where K = 0, ±1 is the FRW curvature index, and m is the mass parameter of the black hole at 
R = (recall that the 5D gravitational potential has R~^ behaviour). The bulk black hole gives 
rise to dark radiation on the brane via its Coulomb effect. The FRW brane moves radially along 
the 5th dimension, with R = a(T), where a is the FRW scale factor, and the junction conditions 
determine the velocity via the Friedmann equation for a |339| 145) . Thus one can interpret the 
expansion of the universe as motion of the brane through the static bulk. In the special case m = 
and da/dT = 0, the brane is fixed and has Minkowski geometry, i.e., the original RS 1-brane 
brane-world is recovered in different coordinates. 

The velocity of the brane is coordinate-dependent, and can be set to zero. We can use Gaussian 
normal coordinates, in which the brane is fixed but the bulk metric is not manifestly static |38j : 



-N^t,y)df+A^{t,y) 



dr^ 



1 - Kr^ 



+ dy\ 



(186) 



Here a{t) = A{t, 0) is the scale factor on the FRW brane at y = 0, and t may be chosen as proper 
time on the brane, so that N{t, 0) = 1. In the case where there is no bulk black hole (m = 0), the 
metric functions are 



N = 



ait) ' 



A = a(t) 



cosh I 



(187) 
(188) 



Again, the junction conditions determine the Friedmann equation. The extrinsic curvature at the 
brane is 

i^^ = diag ("^,4^,4^,4^). (189) 

brane 



N' A' A' A 
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Then, by Equation ([68]) . 



A 



-f{2p + ?,p-\), (190) 

brane O 

^ =-f(p + A). (191) 

brane D 



The field equations yield the first integral [3S] 

{A^f-^A^ + ^A' + ,n = Q, (192) 

where m is constant. Evaluating this at the brane, using Equation (|19ip . gives the modified 
Friedmann equation (|194p . 

The dark radiation carries the imprint on the brane of the bulk gravitational field. Thus we 
expect that Sin, for the Friedmann brane contains bulk metric terms evaluated at the brane. In 
Gaussian normal coordinates (using the field equations to simplify). 



A" 



, S-, = -(\£%]S',. (193) 

brane Z \ o 



Either form of the cosmological metric. Equation (|184p or (|186p . may be used to show that 5D 
gravitational wave signals can take "short-cuts" through the bulk in travelling between points A 
and B on the brane [511 12121 [55] . The travel time for such a graviton signal is less than the time 
taken for a photon signal (which is stuck to the brane) from A to B. 

Instead of using the junction conditions, we can use the covariant 3D Gauss-Codazzi equa- 
tion (|140p to find the modified Friedmann equation: 



on using Equation (|124p . where 



K 



2 



m='—psoa^. (195) 



The covariant Raychauhuri equation p29p yields 



H = -^{p + p)(l + ^)-2^ + ^, (196) 



2 ^'^ V A 

which also follows from differentiating Equation (|194p and using the energy conservation equation. 

When the bulk black hole mass vanishes, the bulk geometry reduces to AdSs, and p£ = 0. 
In order to avoid a naked singularity, we assume that the black hole mass is non-negative, so 
that p£Q > 0. (By Equation p85p . it is possible to avoid a naked singularity with negative m 
when K = —1, provided \m\ < i'^/i.) This additional effective relativistic degree of freedom is 
constrained by nucleosynthesis and CMB observations to be no more than ^ 5% of the radiation 
energy density [1721 Hi, ,209, ,47] : 

— < 0.05. (197) 

Prad nuc 

The other modification to the Hubble rate is via the high-energy correction p/X. In order to 
recover the observational successes of general relativity, the high-energy regime where significant 
deviations occur must take place before nucleosynthesis, i.e., cosmological observations impose the 
lower limit 

A > (1 MeV)^ ^ A/5 > 10** GeV. (198) 
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This is much weaker than the hmit imposed by table-top experiments, Equation (|42[) . Since /\ 
decays as during the radiation era, it will rapidly become negligible after the end of the 
high-energy regime, p = X. 

If P£ = and K ~ = A, then the exact solution of the Friedmann equations for w ~ p/p = 
const, is |38] 



a = const. X + tx = , -^^ < (1 + i(;)-HO-^ s, (199) 

V3^(l + u') ' 

where w > —1. If p£ 7^ (but A' = = A), then the solution for the radiation era (lo = i) is [21] 

a = const. X [t{t -f tx)]"\ = , ■ (200) 

For i ^ we recover from Equations (|199p and (|200p the standard behaviour, a oc ^2/3(111+1)^ 
whereas for t tx, we have the very different behaviour of the high-energy regime, 

/9>A ^ acx (201) 

When w = — 1 we have p = po from the conservation equation. If AT = = A, we recover the 
de Sitter solution for p£ = and an asymptotically de Sitter solution for p£ > 0: 



a = aoexp[i7o(t-to)], ^" = ''V Y + |x) ' /'^ = (202) 

= . ^ sinh[2i/o(t - to)] for ps > 0. (203) 

V Ho 

A qualitative analysis of the Friedmann equations is given in [551 IHZ] • 



5.1 Brane- world inflation 

In 1-brane RS-type brane-worlds, where the bulk has only a vacuum energy, inflation on the brane 
must be driven by a 4D scalar field trapped on the brane. In more general brane-worlds, where the 
bulk contains a 5D scalar field, it is possible that the 5D field induces inflation on the brane via its 
effective projection [232] [TMl ITTn [37T] [TMI ITMl liTOl 1421 IMTI ^ 

More exotic possibilities arise from the interaction between two branes, including possible 
collision, which is mediated by a 5D scalar field and which can induce either inflation [134) 1221) or a 
hot big-bang radiation era, as in the "ekpyrotic" or cyclic scenario [230 1 [2T6 l IMS ] l406 l [275 1 [3T9 1 I4T5 ] . 
or in colliding bubble scenarios [101 11571 1158] . (See also [IHl [Ml 1301] for colliding branes in an 
M theory approach.) Here we discuss the simplest case of a 4D scalar field with potential V{4>) 
(see [289] for a review). 

High-energy brane-world modifications to the dynamics of inflation on the brane have been 
investigated [STOl ^17} M\ EQHI [3221 [Ml [Ml [2871 [HI [35l [36l [3301 UMl [263 ISMl [309] . Essentially, 
the high-energy corrections provide increased Hubble damping, since p ^ X implies that H is 
larger for a given energy than in 4D general relativity. This makes slow-roll inflation possible even 
for potentials that would be too steep in standard cosmology [3101 [99l (Sill (3721 [Mil [288] [206] . 

The field satisfies the Klein-Gordon equation 

4> + 3H(i)+V'{(t)) ^ 0. (204) 

In 4D general relativity, the condition for inflation, d > 0, is (j)^ < V{(t)), i.e., p < —^p, where 
p = ^(j)"^ + V and p = ^(j)'^ — V. The modified Friedmann equation leads to a stronger condition 
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for inflation: Using Equation (|194|) . with m = = A = A', and Equation (|204|) . we find that 



a > 



w < 



"l + 2p/A" 
. 1 + p/A _ 



(205) 



where the square brackets enclose the brane correction to the general relativity result. As p/A — > 0, 
the 4D result w < — ^ is recovered, but for p > X, w must be more negative for inflation. In the 
very high-energy limit p/X — >■ oo, we have w < — |. When the only matter in the universe is a 
self-interacting scalar field, the condition for infiation becomes 

V 



r -v + 



A 



1, 



-V 



<o, 



(206) 



which reduces to (P < V when 



+ F< A. 



In the slow-roll approximation, we get 



iH' 



V_ 
2A 



(207) 
(208) 



The brane- world correction term V/X in Equation (|207p serves to enhance the Hubble rate for 
a given potential energy, relative to general relativity. Thus there is enhanced Hubble 'friction' 
in Equation (|208p . and brane- world effects will reinforce slow- roll at the same potential energy. 
We can see this by defining slow-roll parameters that reduce to the standard parameters in the 
low-energy limit: 

(209) 
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Self-consistency of the slow-roll approximation then requires e, \r]\ <^ 1. At low energies, ^ ^ A, 
the slow-roll parameters reduce to the standard form. However at high energies, V X, the extra 
contribution to the Hubble expansion helps damp the rolling of the scalar field, and the new factors 
in square brackets become « X/V: 



4A 

V 



2A 

v 



(211) 



where egr,??gr are the standard general relativity slow-roll parameters. In particular, this means 
that steep potentials which do not give infiation in general relativity, can inflate the brane-world 
at high energy and then naturally stop inflating when V drops below A. These models can be 
constrained because they typically end inflation in a kinetic-dominated regime and thus generate 
a blue spectrum of gravitational waves, which can disturb nucleosynthesis [99 } 13141 13721 1349. ,288] . 
They also allow for the novel possibility that the inflaton could act as dark matter or quintessence 
at low energies ^ [Hil [372] IM9l IMSl l8l [329] IMH [56] [383 ] . 

The number of e-folds during inflation, N = J Hdt, is, in the slow-roll approximation. 
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(212) 
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Brane- world effects at high energies increase the Hubble rate by a factor y/2A, yielding more 
inflation between any two values of 4> for a given potential. Thus we can obtain a given number of 
e- folds for a smaller initial inflaton value 4>i. For V ^ X, Equation (|212p becomes 



N : 



1287r^ 
3M| 



(213) 



The key test of any modified gravity theory during inflation will be the spectrum of pertur- 
bations produced due to quantum fluctuations of the flelds about their homogeneous background 
values. We will discuss brane- world cosmological perturbations in the next Section [SI In general, 
perturbations on the brane are coupled to bulk metric perturbations, and the problem is very 
complicated. However, on large scales on the brane, the density perturbations decouple from the 
bulk metric perturbations j305l 12731 11781 1149] . For 1-brane RS-type models, there is no scalar 
zero-mode of the bulk graviton, and in the extreme slow-roll (de Sitter) limit, the massive scalar 
modes are heavy and stay in their vacuum state during inflation [149] . Thus it seems a reasonable 
approximation in slow-roll to neglect the KK effects carried by Ef^i, when computing the density 
perturbations. 
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Figure 5: The relation between the inflaton mass m/Mi (M4 = Mp) and the brane tension 
(A/MI)^/'* necessary to satisfy the COBE constraints. The straight line is the approximation 
used in Equation p20|) . which at high energies is in excellent agreement with the exact solution, 
evaluated numerically in slow-roll. (Figure taken from |310) .) 

To quantify the amplitude of scalar (density) perturbations we evaluate the usual gauge- 
invariant quantity 

C = 7^ - -<5p, (214) 
P 

which reduces to the curvature perturbation TZ on uniform density hypersurfaces {5p = 0). This is 
conserved on large scales for purely adiabatic perturbations as a consequence of energy conserva- 
tion (independently of the field equations) |428J . The curvature perturbation on uniform density 
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hypersurfaces is given in terms of the scalar field fluctuations on spatially flat hypersurfaces S(f) by 



(215) 



If one makes the assumption that backreaction due to metric perturbations in the bulk can be 
neglected, the field fluctuations at Hubble crossing (k — aH) in the slow-roll limit are given 
by {Scfy^} sa {H/2'!r) , a result for a massless field in de Sitter space that is also independent 
of the gravity theory j428) . For a single scalar field the perturbations are adiabatic and hence 
the curvature perturbation C, can be related to the density perturbations when modes re-enter the 
Hubble scale during the matter dominated era which is given by = 4(C^)/25. Using the slow-roll 
equations and Equation (|215p . this gives 
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(216) 



k=aH 



Thus the amplitude of scalar perturbations is increased relative to the standard result at a fixed 
value of 4> for a given potential. 

A crucial assumption is that backreaction due to metric perturbations in the bulk can be 
neglected. In the extreme slow-roll limit this is necessarily correct because the coupling between 
inflaton fluctuations and metric perturbations vanishes; however, this is not necessarily the case 
when slow- roll corrections are included in the calculation. Previous work |252l 12551 1257] has shown 
that such bulk effects can be subtle and interesting (see also |110[ 1114] for other approaches). In 
particular, subhorizon inflaton fluctuations on a brane excite an infinite ladder of Kaluza-Klein 
modes of the bulk metric perturbations at first order in slow-roll parameters, and a naive slow-roll 
expansion breaks down in the high-energy regime once one takes into account the backreaction 
of the bulk metric perturbations, as confirmed by direct numerical simulations [201] . However, 
an order-one correction to the behaviour of inflaton fluctuations on subhorizon scales does not 
necessarily imply that the amplitude of the inflaton perturbations receives corrections of order one 
on large scales; one must consistently quantise the coupled brane inflaton fluctuations and bulk 
metric perturbations. This requires a detailed analysis of the coupled brane-bulk system |70 1I254] . 

It was shown that the coupling to bulk metric perturbations cannot be ignored in the equations 
of motion. Indeed, there are order-unity differences between the classical solutions without coupling 
and with slow-roll induced coupling. However, the change in the amplitude of quantum-generated 
perturbations is at next-to-leading order |254] because there is still no mixing at leading order 
between positive and negative frequencies when scales observable today crossed the horizon, so the 
Bogoliubov coefficients receive no corrections at leading order. The amplitude of perturbations 
generated is also subject to the usual slow-roll corrections on super-horizon scales. The next-order 
slow-roll corrections from bulk gravitational perturbations are calculated in [256] and they are the 
same order as the usual Stewart-Lyth correction [407]. These results also show that the ratio of 
tensor-to-scalar perturbation amplitudes are not influenced by brane-bulk interactions at leading 
order in slow-roll. It is remarkable that the predictions from inflation theories should be so robust 
that this result holds in spite of the leading-order change to the solutions of the classical equations 
of motion. 

The scale-dependence of the perturbations is described by the spectral tilt 
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din ^2 
din A: 



-4e + 277, 



(217) 



where the slow-roll parameters are given in Equations ([209]) and (]210j) . Because these slow-roll 
parameters are both suppressed by an extra factor X/V a,i high energies, we see that the spectral 
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index is driven towards the Harrison-Zel'dovich spectrum, Hs — > 1, as V/X — >■ oo; however, as ex- 
plained below, this does not necessarily mean that the branc- world case is closer to scale-invariance 
than the general relativity case. 

As an example, consider the simplest chaotic inflation model V — ^'ni?(jr' . Equation (|212p gives 
the integrated expansion from to as 

The new high-energy term on the right leads to more inflation for a given initial inflaton value (^j. 

The standard chaotic inflation scenario requires an inflaton mass m ~ 10^^ GeV to match the 
observed level of anisotropics in the cosmic microwave background (see below). This corresponds 
to an energy scale ~ 10^^ GeV when the relevant scales left the Hubble scale during inflation, 
and also to an inflaton field value of order 3Mp. Chaotic inflation has been criticised for requiring 
super-Planckian field values, since these can lead to nonlinear quantum corrections in the potential. 

If the brane tension A is much below 10^^ GeV, corresponding to M5 < 10^^ GeV, then the 
terms quadratic in the energy density dominate the modified Friedmann equation. In particular 
the condition for the end of inflation given in Equation (|206p becomes (f^ < ^V. In the slow-roll 
approximation (using Equations (|207p and (|208p ) (j) ~ ^^^l/STrt/f, and this yields 

In order to estimate the value of when scales corresponding to large-angle anisotropics on the 
microwave background sky left the Hubble scale during inflation, we take A^cobe ~ 55 in Equa- 
tion (|218p and (pi = (f>cnd- The second term on the right of Equation (|218p dominates, and we 
obtain 

Imposing the COBE normalization on the curvature perturbations given by Equation (j216p requires 

Substituting in the value of ^cobo given by Equation (|220p shows that in the limit of strong brane 
corrections, observations require 

m « 5 X 10~^ M5, 0cobc « 3 X 10^ M5. (222) 

Thus for A/,5 < 10^^ GeV, chaotic inflation can occur for field values below the 4D Planck scale, 
0cobc < Mp, although still above the 5D scale M^. The relation determined by COBE constraints 
for arbitrary branc tension is shown in Figure [SJ together with the high-energy approximation used 
above, which provides an excellent fit at low brane tension relative to M4. 

It must be emphasized that in comparing the high-energy brane-world case to the standard 
4D case, we implicitly require the same potential energy. However, precisely because of the high- 
energy effects, large-scale perturbations will be generated at different values of V than in the 
standard case, specifically at lower values of V, closer to the reheating minimum. Thus there are 
two competing effects, and it turns out that the shape of the potential determines which is the 
dominant effect j286) . For the quadratic potential, the lower location on V dominates, and the 
spectral tilt is slightly further from scale invariance than in the standard case. The same holds for 
the quartic potential. Data from WMAP and 2dF can be used to constrain inflationary models via 
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Figure 6: Constraints from WMAP data on inflation models with quadratic and quartic potentials, 
where R is the ratio of tensor to scalar amplitudes and n is the scalar spectral index. The high 
energy (H.E.) and low energy (L.E.) limits are shown, with intermediate energies in between, and 
the l-cr and 2-cr contours are also shown. (Figure taken from |286) .) 
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their deviation from scale invariance, and the high-energy brane- world versions of the quadratic and 
quartic potentials are thus under more pressure from data than their standard counterparts [286] , 
as shown in Figure [S] 

Other perturbation modes have also been investigated: 

• High-energy inflation on the brane also generates a zero-mode (4D graviton mode) of tensor 
perturbations, and stretches it to super-Hubble scales, as will be discussed below. This zero- 
mode has the same qualitative features as in general relativity, remaining frozen at constant 
amplitude while beyond the Hubble horizon. Its amplitude is enhanced at high energies, 
although the enhancement is much less than for scalar perturbations [274] : 
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Equation (|224p means that brane-world effects suppress the large-scale tensor contribution 
to CMB anisotropies. The tensor spectral index at high energy has a smaller magnitude than 
in general relativity, 

nt ^ -3e, (225) 
but remarkably the same consistency relation as in general relativity holds |206| : 

nt = -2^. (226) 

This consistency relation persists when Zi symmetry is dropped [207] (and in a two-brane 
model with stabilized radion |173| ). It holds only to lowest order in slow-roll, as in general 
relativity, but the reason for this |384] and the nature of the corrections [64] are not settled. 

The massive KK modes of tensor perturbations remain in the vacuum state during slow-roll 
inflation [2741 1177) . The evolution of the supcr-Hubblc zero mode is the same as in general 
relativity, so that high-energy brane-world effects in the early universe serve only to rescale 
the amplitude. However, when the zero mode re-enters the Hubble horizon, massive KK 
modes can be excited. 

• Vector perturbations in the bulk metric can support vector metric perturbations on the brane, 
even in the absence of matter perturbations (see the next Section IS]). However, there is no 
normalizablc zero mode, and the massive KK modes stay in the vacuum state during brane- 
world inflation |52| . Therefore, as in general relativity, we can neglect vector perturbations 
in inflationary cosmology. 

Brane-world effects on large-scale isocurvature perturbations in 2-field inflation have also been 
considered [17]. Brane-world (p)reheating after inflation is discussed in [4171 1432] [9] I418| l96]. 



5.2 Brane-world instanton 

The creation of an inflating brane-world can be modelled as a de Sitter instanton in a way that 
closely follows the 4D instanton, as shown in |155| . The instanton consists of two identical patches of 
AdSs joined together along a de Sitter brane (dS4) with compact spatial sections. The instanton 
describes the "birth" of both the inflating brane and the bulk spacetime, which are together 
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"created from nothing", i.e., the point at the south pole of the de Sitter 4-sphere. The Euchdean 
AdSs metric is 



dx + sin X 



(3) 



(227) 



where ^^^(3) is a 3-sphere, and r < vq. The Euchdean instanton interpolates between r = 
("nothing") and r = tq (the created universe), which is a spherical brane of radius 

= £smh{ro/e). (228) 

After creation, the brane-world evolves according to the Lorentzian continuation, x ^ iHot + tt/2, 



-dr 



HQ^cosh^{Hot) dn 



(3) 



(see Figure [7]). 



r = 




(229) 



Figure 7: Brane-world instanton. (Figure taken from |155] .) 



5.3 Models with non-empty bulk 

The single-brane cosmological model can be generalized to include stresses other than A5 in the 
bulk: 

• The simplest example arises from considering a charged bulk black hole, leading to the 
Reissner-Nordstrom AdSs bulk metric [22]. This has the form of Equation (jl84p . with 



(230) 



where q is the "electric" charge parameter of the bulk black hole. The metric is a solution of 
the 5D Einstein-Maxwell equations, so that '^^^Tab in Equation (|50p is the energy- momentum 
tensor of a radial static 5D "electric" field. In order for the field lines to terminate on 
the boundary brane, the brane should carry a charge ~q. Since the RNAdSs metric is 4- 
isotropic, it is still possible to embed a FRW brane in it, which is moving in the coordinates 
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of Equation (|184|) . The effect of the black hole charge on the brane arises via the junction 
conditions and leads to the modified Friedmann equation |22j . 



^ =yKi + l^) + ^-!a + 3^-:;.- (231) 



The field lines that terminate on the brane imprint on the brane an effective negative energy 
density —8(7^/(^^0®), which redshifts like stiff matter [w = 1). The negativity of this term 
introduces the possibility that at high energies it can bring the expansion rate to zero and 
cause a turn-around or bounce (but see [205] for problems with such bounces). 

Apart from negativity, the key difference between this "dark stiff matter" and the dark 
radiation term m/a^ is that the latter arises from the bulk Weyl curvature via the 
tensor, while the former arises from non-vacuum stresses in the bulk via the J^^^, tensor in 
Equation (|69p . The dark stiff matter does not arise from massive KK modes of the graviton. 



• Another example is provided by the Vaidya-AdS5 metric, which can be written after trans- 
forming to a new coordinate u = T + / dR/F in Equation (|184p . so that v = const, are null 
surfaces, and 

(^ks^ = ~F{R, v)dv^ + 2dv dR + i?2 Y^^^ + r^dn'^^ , (232) 

m.) = A- + ^-^f>. (233) 

This model has a moving FRW brane in a 4-isotropic bulk (which is not static), with either a 
radiating bulk black hole (dm/dv < 0), or a radiating brane {dm/dv > 0) [771 12801 [?7^ 1282] . 
The metric satisfies the 5D field equations (|50|) with a null-radiation energy-momentum 
tensor, 

(■5)rAi3 = VfcAfciJ, kAk^ = 0, kAU^ = 1, (234) 

where ip oc dm/dv. It follows that 

J^Hf = K^'^-tphfj.iy (235) 

In this case, the same effect, i.e., a varying mass parameter m, contributes to both 5^^^, and 
J^fj^i, in the brane field equations. The modified Friedmann equation has the standard 1-brane 
RS-type form, but with a dark radiation term that no longer behaves strictly like radiation: 



K 



2 



^^^yHl + i^J+^+3^-S- (236) 
By Equations ^ and ((234| . we arrive at the matter conservation equations, 

VT^, = ~2^u^. (237) 

This shows how the brane loses {tp > 0) or gains {ip < 0) energy in exchange with the bulk 
black hole. For an FRW brane, this equation reduces to 

p + 3H{p + p) ^ ^2i;. (238) 

The evolution of m is governed by the 4D contracted Bianchi identity, using Equation (|235p : 



\7t^£ — 



(V^-f 61/') -3k2i^] u^ + ^K^V^^. (239) 
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For an FRW brane, this yields 
where = 3m(t)/(K^a'*). 



(240) 



• A more comphcated bulk metric arises when there is a self-interacting scalar field $ in the 
bulk [313 [2IJ [321 [Till [2761 [1451 SI]. In the simplest case, when there is no coupling between 
the bulk field and branc matter, this gives 



(5; 



Tab = $,A$,i3 



(5), 



gAB 



where $(a;,?/) satisfies the 5D Klein-Gordon equation, 

(■5)n$_y'($) ^0. 
The junction conditions on the field imply that 

dy^{x,Q) = 0. 



(241) 



(242) 



(243) 



Then Equations (j74p and (|24ip show that matter conservation continues to hold on the brane 
in this simple case: 

V^T^, = 0. (244) 

From Equation p4ip one finds that 



40.^0,. - g^u (^3y(</)) + |g"''0,a0,;j 



where 



(/.(x) = $(x,0), 
so that the modified Friedmann equation becomes 
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-6' + V{4>) 



(245) 
(246) 

(247) 



When there is coupling between brane matter and the bulk scalar field, then the Friedmann 
and conservation equations are more complicated |313l [^Tl 13241 11441 12761 11451 [3H] • 
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6 Brane- World Cosmology: Perturbations 

The background dynamics of brane-world cosmology are simple because the FRW symmetries 
simplify the bulk and rule out nonlocal effects. But perturbations on the brane immediately 
release the nonlocal KK modes. Then the 5D bulk perturbation equations must be solved in order 
to solve for perturbations on the brane. These 5D equations are partial differential equations for 
the 3D Fourier modes, with both initial and boundary conditions needed. 

The theory of gauge-invariant perturbations in brane-world cosmology has been extensively 

investigated and developed ^ 'm^ m \m[ m \m[ \m\ \m\ \m\ \TW j[7T 

um nm m7\ wm m]\ imsi [420. 2321 \m\ itmi um \m\ m nm ms. .388. .371. ,2571 m 

M nm mi m MH ESI MB ISMI mm and is qualitatively well understood. The key task is 
integration of the coupled brane-bulk perturbation equations. Special cases have been solved, 
where these equations effectively decouple |273l [MJ 12841 1283) , and approximation schemes have 
been developed [inilll3Il[Mnilin3lini[M71IMl[SIl[ini[I33 more general 

cases where the coupled system must be solved. Below we will also present the results of full 
numerical integration of the 5D perturbation equations in the RS case. 

From the brane viewpoint, the bulk effects, i.e., the high-energy corrections and the KK modes, 
act as source terms for the brane perturbation equations. At the same time, perturbations of mat- 
ter on the brane can generate KK modes (i.e., emit 5D gravitons into the bulk) which propagate 
in the bulk and can subsequently interact with the brane. This nonlocal interaction amongst the 
perturbations is at the core of the complexity of the problem. It can be elegantly expressed via 
integro-differential equations |334l 1336] . which take the form (assuming no incoming 5D gravita- 
tional waves) 

Ak{t) = J dt'g{t,t')Bk{t'), (248) 

where Q is the bulk retarded Green's function evaluated on the brane, and Ak,Bk are made 
up of brane metric and matter perturbations and their (brane) derivatives, and include high- 
energy corrections to the background dynamics. Solving for the bulk Green's function, which then 
determines G, is the core of the 5D problem. 

We can isolate the KK anisotropic stress tt^^ as the term that must be determined from 5D 
equations. Once tt^^, is determined in this way, the perturbation equations on the brane form a 
closed system. The solution will be of the form (expressed in Fourier modes): 

^f(Ocx j dt' g{t,t')Fk{t'), (249) 

where the functional will be determined by the covariant brane perturbation quantities and their 
derivatives. It is known in the case of a Minkowski background [377], but not in the cosmological 
case. 

The KK terms act as source terms modifying the standard general relativity perturbation 
equations, together with the high-energy corrections. For example, the linearization of the shear 
propagation equation (jl3ip yields 

2 2 2 

a^,^ + 2Ha^,^ + E^,^ - yTT^^ - V^^.A^) = yTr^^ " + 3w)^7r^^. (250) 

In 4D general relativity, the right hand side is zero. In the brane-world, the first source term on the 
right is the KK term, and the second term is the high-energy modification. The other modification 
is a straightforward high-energy correction of the background quantities H and p via the modified 
Friedmann equations. 

As in 4D general relativity, there are various different, but essentially equivalent, ways to 
formulate linear cosmological perturbation theory. First we describe the covariant brane-based 
approach. 
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6.1 1 + 3-covariant perturbation equations on the brane 

In the 1 + 3-covariant approach j305l 12841 1307] . perturbative quantities are projected vectors, 
= V^^) , and projected symmetric tracefree tensors, W^i/ = Wi^^^^^ , which are gauge-invariant 
since they vanish in the background. These are decomposed into (3D) scalar, vector, and tensor 
modes as 

V^,=V,,V + V^, (251) 
W^, - V(^V,)Ty + V^^iy.) -t- W^,, (252) 

where W^i/ = W(^fj_^^ and an ovcrbar denotes a (3D) transverse quantity, 

V^V;, =0 = V'^VK^,. (253) 

In a local inertial frame comoving with u^, i.e., — (1,0), all time components may be set to 
zero: = (0, V,), Wo^ = 0, = (0, V,). 

Purely scalar perturbations are characterized by the fact that vectors and tensors are derived 
from scalar potentials, i.e., 

= W,, = W^, = 0. (254) 

Scalar perturbative quantities are formed from the potentials via the (3D) Laplacian, e.g., V = 
V^V^T^ = V^l^. Purely vector perturbations are characterized by 

V^ = V^, W^M- - V(^VF,), curlV^/ = -2/w^, (255) 

where uj^ is the vorticity, and purely tensor by 

V^./ = = T/^, W^^. = W^p.. (256) 

The KK energy density produces a scalar mode V ^ps (which is present even if = in the 
background). The KK momentum density carries scalar and vector modes, and the KK anisotropic 
stress carries scalar, vector, and tensor modes: 

+ (257) 
T^iu = V(^V,)^^ + V(^^f^ -f Til,. (258) 

Linearizing the conservation equations for a single adiabatic fluid, and the nonlocal conservation 
equations, we obtain 

p + e(p + p) = 0, (259) 
c'V^ + (p + p)Ap = 0, (260) 

h + \QP£ + = 0, (261) 



ql + mql + iv^p£ + ^p,A^ + V'^^^, = -I^l^V^p. (262) 
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Linearizing the remaining propagation and constraint equations leads to 
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(263) 
(264) 
(265) 

(266) 

(267) 
(268) 
(269) 
(270) 
(271) 

(272) 



Equations p59p . p6ip . and ([263]) do not provide gauge-invariant equations for perturbed quanti- 
ties, but their spatial gradients do. 

These equations are the basis for a 1 -I- 3-covariant analysis of cosmological perturbations from 
the brane observer's viewpoint, following the approach developed in 4D general relativity (for a 
review, see |138j ). The equations contain scalar, vector, and tensor modes, which can be separated 
out if desired. They are not a closed system of equations until tt^^^ is determined by a 5D analysis 
of the bulk perturbations. An extension of the 1 -I- 3-covariant perturbation formalism to 1 -I- 4 
dimensions would require a decomposition of the 5D geometric quantities along a timelike extension 

into the bulk of the brane 4- velocity field m'^, and this remains to be done. The 1 -I- 3-covariant 
perturbation formalism is incomplete until such a 5D extension is performed. The metric-based 
approach does not have this drawback. 



6.2 Metric-based perturbations 

An alternative approach to branc-world cosmological perturbations is an extension of the 4D 
metric-based gauge- invariant theory |2371 1331] . A review of this approach is given in |53[ 1365] . In 
an arbitrary gauge, and for a flat FRW background, the perturbed metric has the form 









Si) 


Na 






{27^% + 2d.,djC - 


V2d^,F^)+hj] 


A\^,P-X^) 




Na 


A^{d,p- 




2v 



where the backgroimd metric functions A,N are given by Equations (|187l I188p . The scalars 
ip,TZ,C,a, represent scalar perturbations. The vectors Si, Fi, and Xi ^^e transverse, so that 
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they represent 3D vector perturbations, and the tensor fij is transverse traceless, representing 3D 
tensor perturbations. 

In the Gaussian normal gauge, the brane coordinate-position remains fixed under perturbation. 



'"rfs' = [gl^'Jix, y) + y)\ dx'^dx'^ + dy\ 

where g^^} is the background metric, Equation (jl86p . In this gauge, we have 

« = = = Xj=0. 
In the 5D longitudinal gauge, one gets 

-B + C = {] = -P + C'. 



(274) 



(275) 



(276) 



In this gauge, and for an AdSs background, the metric perturbation quantities can all be expressed 
in terms of a "master variable" VI which obeys a wave equation |334l 1336] . In the case of scalar 
perturbations, we have for example 



7^ = 



(277) 



with similar expressions for the other quantities. All of the metric perturbation quantities are 
determined once a solution is found for the wave equation 



1 



n 



P \ N 



6 v42 / A3 



N 
A^ 



n' 



(278) 



The junction conditions (|68p relate the ofF-brane derivatives of metric perturbations to the 
matter perturbations: 



dy 5g^^ 



where 



<5r% - -Sp, 

For scalar perturbations in the Gaussian normal gauge, this gives 



dy^{x,Q) = -^{25p + ?,5p), 
dyB{x, 0) = nlSp, 



Sir, 



dyC{x,0) 
dyTZ{x,0) 

where Sn is the scalar potential for the matter anisotropic stress. 



-f5p~d,d'C{x,0), 



(279) 



(280) 
(281) 
(282) 



(283) 
(284) 

(285) 
(286) 

(287) 
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The perturbed KK energy-momentum tensor on the brane is given by 

(5f % = i^'^hs, 



5E^ 



(288) 
(289) 

(290) 



The evolution of the bulk metric perturbations is determined by the perturbed 5D field equa- 
tions in the vacuum bulk, 



8 (5)G" 



0. 



(291) 



Then the matter perturbations on the brane enter via the perturbed junction conditions (j279p . 
For example, for scalar perturbations in Gaussian normal gauge, we have 



V'-27^' 



2iV2 



.A N 

'a^n 



B' 



(292) 



For tensor perturbations (in any gauge), the only nonzero components of the perturbed Einstein 
tensor are 



(293) 



In the following, we will discuss various perturbation problems, using either a 1 -|- 3-covariant 
or a metric-based approach. 



6.3 Density perturbations on large scales 

In the covariant approach, we define matter density and expansion (velocity) perturbation scalars, 
as in 4D general relativity. 



z = a^v^e. 



Then we can define dimensionless KK perturbation scalars |305) . 



U 



Q = -vV, 

P 



(294) 



(295) 



where the scalar potentials and tt^ are defined by Equations P571I258|) . The KK energy density 
(dark radiation) produces a scalar fiuctuation U which is present even ii pg = in the background, 
and which leads to a non-adiabatic (or isocurvature) mode, even when the matter perturbations 
arc assumed adiabatic [178) . We define the total effective dimensionless entropy S'tot via 



2t~72 2 2t~72 

a V ptot - Ctotfl V ptot. 



Ptot S'tot 

where c^q^ = ptot/ptot is given in Equation (|109p . Then 
e _ 9[c2-i + (|+^ + c2)p/A] 



[3(1 + w){l + p/\) + Ap£/p][iw + 3(1 -f 2w)p/2A + ps/p] 



3 p 



(296) 



(297) 



If P£ = in the background, then U is an isocurvature mode: S'tot oc {\+w)U . This isocurvature 
mode is suppressed during slow-roll inflation, when 1 + « 0. 
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If P£ 7^ in the background, then the weighted difference between U and A determines the 
isocurvaturc mode: Stot oc (4p£/3p)A — (1 + w)U . At very high energies, p ^ X, the entropy is 
suppressed by the factor X/p. 

The density perturbation equations on the brane are derived by taking the spatial gradients of 
Equations (fSMj) . (f26T|) . and (|2H3l) . and using Equations ([260]) and (|262|) . This leads to [178] 



A = 3wHA - (1 + 
1 



1 + (4 + 3w) 



id 



1+w 



P£_ 
P 



U =(3w- 1)HU 



1 



If 



P 



l)HQ-^U-laV^n+^ 
6a 6 ip, 



4c2 



3 1 + UI - 

P >^ 



(298) 
(299) 

(300) 

(301) 



The KK anisotropic stress term 11 occurs only via its Laplacian, V^II. If we can neglect this 
term on large scales, then the system of density perturbation equations closes on super-Hubble 
scales |305j . An equivalent statement applies to the large-scale curvature perturbations |273) . KK 
effects then introduce two new isocurvaturc modes on large scales (associated with U and Q), and 
they modify the evolution of the adiabatic modes as well |178[ 1284) . 

Thus on large scales the system of brane equations is closed, and we can determine the density 
perturbations without solving for the bulk metric perturbations. 

We can simplify the system as follows. The 3-Ricci tensor defined in Equation (|140p leads to a 
scalar covariant curvature perturbation variable. 



C = a'^V^R^ = -Aa^HZ + 2K^a^p (i + -^\a + 2K^a^pU. 

V 2A / 

It follows that C is locally conserved (along flow lines): 

C = Co, Co = 0. 
We can further simplify the system of equations via the variable 

$ = K^a^pA. 



(302) 



(303) 



(304) 



This should not be confused with the Bardeen metric perturbation variable , although it is the 
covariant analogue of in the general relativity limit. In the brane- world, high-energy and KK 
effects mean that is a complicated generalization of this expression |284) involving 11, but the 
simple $ above is still useful to simplify the system of equations. Using these new variables, we 
find the closed system for large-scale perturbations: 
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H 
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3a^Hp 
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If there is no dark radiation in the background, p£ ^ 0, then 



U ~ Uq exp 



(1 - 3w)dN 



(305) 
(306) 

(307) 
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Figure 8: The evolution of the covariant variable $, defined in Equation (|304p (and not to be 
confused with the Bardeen potential) , along a fundamental world-line. This is a mode that is well 
beyond the Hubble horizon at iV = 0, about 50 e-folds before inflation ends, and remains super- 
Hubblc through the radiation era. A smooth transition from inflation to radiation is modelled by 
= i [(2 — |e) tanh(iV — 50) — (1 — |e)] , where e is a smaU positive parameter (chosen as e = 0.1 in 
the plot). Labels on the curves indicate the value of po/X, so that the general relativistic solution 
is the dashed curve {po/X = 0). (Figure taken from |178] .) 
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and the above system reduces to a single equation for $. At low energies, and for constant w, the 
non-decaying attractor is the general relativity solution, 



3(1 + 



At very high energies, for w > — i, we get 



2(5 + 3w) 



(308) 



3 A . 

$high^ T^ — (1 + w) 

^ Po 



Co 



2[/o 



7 + Gw 5 + 



(309) 



where Uq = K^agPoUo, so that the isocurvature mode has an influence on $. Initially, $ is 
suppressed by the factor X/ po, but then it grows, eventually reaching the attractor value in Equa- 
tion (|308p . For slow-roll inflation, when 1 + w ^ e, with < e ^ 1 and i7~^|e| ~ |e'| ^ 1, we get 



high 



3 A 
-e— Coe 
2 Po 



3eN 



(310) 



where A'' = ln(a/ao), so that $ has a growing-mode in the early universe. This is different from 
general relativity, where $ is constant during slow-roll inflation. Thus more amplification of $ can 
be achieved than in general relativity, as discussed above. This is illustrated for a toy model of 
inflation-to-radiation in Figure [8] The early (growing) and late time (constant) attractor solutions 
are seen explicitly in the plots. 

The presence of dark radiation in the background introduces new features. In the radiation era 
{w = i); the non-decaying low-energy attractor becomes [188] 



Co 



(1-a), 



a = — < 0.05. 



(311) 
(312) 



The dark radiation serves to reduce the final value of $, leaving an imprint on $, unlike the p£ 
case. Equation p08|) . In the very high energy limit. 
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(313) 



Thus $ is initially suppressed, then begins to grow, as in the no-dark-radiation case, eventually 
reaching an attractor which is less than the no-dark-radiation attractor. This is confirmed by the 
numerical integration shown in Figure |9l 



6.4 Curvature perturbations and the Sachs— Wolfe effect 



The curvature perturbation TZ on uniform density surfaces is defined in Equation (j273p . The 
associated gauge-invariant quantity 

may be deflned for matter on the brane. Similarly, for the Weyl "fluid" if 7^ in the background, 
the curvature perturbation on hypersurfaces of uniform dark energy density is 



4p£ 



(315) 



53 




Figure 9: The evolution of $ in the radiation era, with dark radiation present in the background. 
(Figure taken from |188| .) 



54 



On large scales, the perturbed dark energy conservation equation is [273] 



(316) 



which leads to 

Cs = 0. (317) 
For adiabatic matter perturbations, by the perturbed matter energy conservation equation, 

{Spy +3H(Sp + Sp) + 3{p + p)'lZ = 0, (318) 



we find 



C = o. 



(319) 



This is independent of brane- world modifications to the field equations, since it depends on en- 
ergy conservation only. For the total, effective fluid, the curvature perturbation is defined as 
follows |273| : li p£ ^ in the background, we have 



Ctot — c + 



4yO£ 



3{p + p){l + p/X)+Ap£ 



iCe ~ C) , 



and ii p£ = in the background, wc get 

Ctot = c + 
Sps - -r, 



Spe 



3{p + p){l + p/X) 



(320) 

(321) 
(322) 



where SC'£ is a constant. It follows that the curvature perturbations on large scales, like the density 
perturbations, can be found on the brane without solving for the bulk metric perturbations. 

Note that Ctot 7^ even for adiabatic matter perturbations; for example, ii p£ = in the 
background, then 

2 1 ^ ^P£ 



Ct 



H 



{p + p){l + p/X)- 



(323) 



The KK effects on the brane contribute a non-adiabatic mode, although Ctot —J' at low energies. 

Although the density and curvature perturbations can be found on super-Hubble scales, the 
Sachs- Wolfe effect requires tt^^ in order to translate from density/curvature to metric perturba- 
tions. In the 4D longitudinal gauge of the metric perturbation formalism, the gauge-invariant 
curvature and metric perturbations on large scales are related by 



Ctot^7e-|(^|-V' 

7^ + V' = -K^a^dTT£, 



(324) 
(325) 



where the radiation anisotropic stress on large scales is neglected, as in general relativity, and Stt£ is 
the scalar potential for tt^^,, equivalent to the covariant quantity II deflned in Equation ()295|) . In 4D 
general relativity, the right hand side of Equation p25p is zero. The (non-integrated) Sachs-Wolfe 
formula has the same form as in general relativity: 



ST 
~T 



= (Crad+^-7^)| 



dec ■ 



(326) 
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The brane-world corrections to the general relativistic Sachs- Wolfe effect are then given by |273j 



where Ss is the KK entropy perturbation (determined by Spg). The KK term Sirg cannot be 
determined by the 4D brane equations, so that ST/T cannot be evaluated on large scales without 
solving the 5D equations. (Equation p27p has been generalized to a 2-branc model, in which the 
radion makes a contribution to the Sachs- Wolfe effect j246j .) 

The presence of the KK (Weyl, dark) component has essentially two possible effects: 

• A contribution from the KK entropy perturbation that is similar to an extra isocurvature 
contribution. 

• The KK anisotropic stress Sirs also contributes to the CMB anisotropics. In the absence 
of anisotropic stresses, the curvature perturbation Ctot would be sufficient to determine the 
metric perturbation TZ and hence the large-angle CMB anisotropies via Equations p24[ 13251 
I326p . However, bulk gravitons generate anisotropic stresses which, although they do not 
affect the large-scale curvature perturbation Ctot, can affect the relation between ^tot: T^-, a-nd 
■0, and hence can affect the CMB anisotropies at large angles. 

A simple phenomenological approximation to Sng on large scales is discussed in |24j . and the 
Sachs- Wolfe effect is estimated as 



ST ^ 






'ln(iin/i4)" 
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^ ^dcc / 


_ln(tcq/t4). 



(328) 



where is the 4D Planck time, and tin is the time when the KK anisotropic stress is induced on 
the brane, which is expected to be of the order of the 5D Planck time. 

A self-consistent approximation is developed in [247] . using the low-energy 2-brane approxi- 
mation gnu USB [MOl HOI 1103] to find an effective 4D form for S^^ and hence for Sirs. This is 
discussed below. In a single brane model in the AdS bulk, full numerical simulations were done to 
find the behaviour of Stts 1691 , as will be discussed in the next subsection. 



6.5 Full numerical solutions 



In order to study scalar perturbations fully, we need to numerically solve the coupled bulk and 
brane equations for the master variable fl. A thorough analysis was done in |69j . For this purpose, 
it is convenient to use the static coordinate where the bulk equation is simple and consider a 
moving FRW brane; 

(5)ds2 ^ r^^^rfj-M^;^^ 4. dz^^^ (329) 



The bulk master variable satisfies the following wave equation (see Equation (|278p ) 







d^n 3 dn 

dz^ z dz 



n. 



(330) 



From the junction condition, O satisfies a boundary condition on the brane 

1 iO\ „ 6pa^ 



1 



0. 



(331) 
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where dn is the derivative orthogonal to the brane 

d 



dn 



OT 



(332) 



and A is the density perturbation in the comoving gauge. The subscript b imphes that the 
quantities are evaluated on the brane. On a brane, A satisfies a wave equation 



drj'^ arj 

A = 6c? - 1 - 8w + 3w^ 
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A = 



k^{i + w)nt 



B = 'ici-Qw~ 4, 



(333a) 
(333b) 



where we consider a perfect fluid on a brane with an equation state w and Cs is a sound speed 
for perturbations. The above ordinal differential equation, the bulk wave equation p30[) and the 
boundary condition p3ip comprise a closed set of equations for A and fie- 
On a brane, we take the longitudinal gauge 



-(1 + 2^lj)df + (1 + 2TZ)S,jdx'dx^ . 



(334) 



Using the expressions for metric perturbations in terms of the master variable Q (see Eq. ([277])) 
and the junction condition Eq. p3ip . ip and TZ are written in terms of A and as 
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2£a3 dry^ ' 
(335b) 



Other quantities of interest are the curvature perturbation on uniform density slices. 



HaV 
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3(1 



ipa^iwX - X- p) 
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Ha dA 



1 + w k^{l + w) drj Ua^ 



f^b, (336) 



where the velocity perturbation V is also written by A and 51b- There are two independent 
numerical codes that can be used to solve for A and fib- The first is the pseudo-spectral (PS) 
method used in |201] and the second is the characteristic integration (CI) algorithm developed in 

Figure [TU] shows the output of the PS and CI codes for a typical simulation of a mode with 
p/X = 50 at the horizon re-enter. As expected we have excellent agreement between the two codes, 
despite the fact that they use different initial conditions. Note that for all simulations, we recover 
that A and C are phase-locked plane waves, 



A(7,)cxcos^, A(7,) « 4C(r/), 



(337) 



at sufficiently late times kr] ^ 1, which is actually the same behaviour as seen in GR. Figure [TT] 
illustrates how the ordinary superhorizon behaviour of perturbations in GR is recovered for modes 
entering the Hubble horizon in the low energy era. We see how A, ip and TZ smoothly interpolate 
between the non-standard high-energy behaviour to the usual expectations in GR. Also shown 
in this plot is the behaviour of the KK anisotropic stress, which steadily decays throughout the 
simulation. These results confirm that at low energies, we recover GR solutions smoothly. 

At high energies p > X, there are two separate effects to consider: First, there is the modifica- 
tion of the universe's expansion at high energies and the 0{p/X) corrections to the perturbative 
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brane scale factor log(a/a,) 

Figure 10: Comparison between typical results of the PS and CI codes for various brane quantities 
(left); and the typical behaviour of the bulk master variable (right) as calculated by the CI method. 
Very good agreement between the two different numerical schemes is seen in the left panel, despite 
the fact that they use different initial conditions. Also note that on subhorizon scales, A and 
^ undergo simple harmonic oscillations, which is consistent with the behaviour in GR. The bulk 
profile demonstrates our choice of initial conditions: We see that the bulk master variable is 
essentially zero during the early stages of the simulation, and only becomes "large" when the mode 
crosses the horizon. Figure taken from |69) . 



equations of motion. Second, there is the effect of the bulk degrees of freedom encapsulated by 
the bulk master variable ft (or, equivalently, the KK fluid To separate out the two effects, 

it is useful to introduce the 4-dimensional effective theory where all 0{p/\) corrections to GR 
are retained, but the bulk effects are removed by artificially setting = 0. In the case of radia- 
tion domination, we obtain equations for the effective theory density contrast A,„f and curvature 
perturbation C,^^f from Equations p36p and p33p with ilb = 0: 

- U 2Xk^ + 4AW ) + ^^^^^^ 

These give a closed set of ODEs on the brane that describe all of the 0{p/X) corrections to GR. 

Since in any given model we expect the primordial value of the curvature perturbation to be 
fixed by inflation, it makes physical sense to normalize the waveforms from each theory such that 
CsD ~ Cet ~ Cgr ~ 1 for a <^ a^,. We can define a set of "enhancement factors" , which are functions 
of k that describe the relative amplitudes of A after horizon crossing in the various theories. 
Let the final amplitudes of the density perturbation with wavenumber k be Csoik), Cret(fc) and 
Cgr(^) for the 5-dimensional, effective and GR theories, respectively, given that the normalization 
CsD ~ Cot ~ Cgr ~ 1 holds. Then, we define enhancement factors as 

^^"^-Wr e-<^'^lf ■ 

It follows that Q„t(fc) represents the 0{p/X) enhancement to the density perturbation, Qs{k) gives 
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Figure 11: The simulated behaviour of a mode on superhorizon scales. On the left we show 
how the A gauge invariant switches from the high-energy behaviour predicted to the familiar GR 
result as the universe expands through the critical epoch. We also show how the KK anisotropic 
stress k\5its steadily decays throughout the simulation, which is typical of all the cases we have 
investigated. On the right, we show the metric perturbations ?/;(= and TZ{= $) as well as the 
curvature perturbation C^. Again, note how the GR result $ « —5' w — 2(/3 is recovered at low 
energy. Figure taken from |69j . 
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the magnification due to KK modes, while Qsoik) gives the total 5-dimensional amplification over 
the GR case. They all increase as the scale is decreased, and that they all approach unity for fc — )■ 0. 
Since Q — 1 implies no enhancement of the density perturbations over the standard result, this 
means we recover general relativity on large scales. For all wavenumbers we see Q„.f > Qs > 1, 
which implies that the amplitude magnification due to the 0{p/\) corrections is always larger 
than that due to the KK modes. Interestingly, the Q-factors appear to approach asymptotically 
constant values for large k: 

S„.f(fc)«3.0, Q£(fc)«2.4, Q5D(fc)«7.1, fc>fcc, (340) 

where kc is the comoving wavenumber of the mode that enters the horizon when H = £~^. 

In cosmological perturbation theory, transfer functions are very important quantities. They 
allow one to transform the primordial spectrum of some quantity set during inflation into the 
spectrum of another quantity at a later time. In this sense, they are essentially the Fourier 
transform of the retarded Green's function for cosmological perturbations. There are many different 
transfer functions one can define, but for our case it is useful to consider a function T(fc) that will 
tell us how the initial spectrum of curvature perturbations T'^"* maps onto the spectrum of density 
perturbations Va at some low energy epoch within the radiation era. It is customary to normalize 
transfer functions such that T(fc; 77) — >■ l(fc — 0), which leads us to the following definition 



-2 



(341) 



Here, Cfc"' is the primordial value of the curvature perturbation and A.k{ri) is the maximum ampli- 
tude of the density perturbation in the epoch of interest. As demonstrated in Figure [TTl we know 
that we recover the GR result in the extreme small scale limit (fc -> 0), which gives the transfer 
function the correct normalization. In the righthand panel of Figure 1121 we show the transfer 
functions derived from GR, the effective theory and the 5-dimensional simulations. As expected, 
the T{k;ri) for each formulation match one another on subcritical scales k < kc- However, on 
supercritical scales we have T^-o > ^lot > Tgr. 

Note that if we are interested in the transfer function at some arbitrary epoch in the low-energy 
radiation regime Ha ^ k^, it is approximately given in terms of the enhancement factor as follows: 

Jl, k<3Ha, 
'''^"^ \(3i7a/fc)2Q,„(fc), k>3Ha, ^^^^^ 

Now, the spectrum of density fluctuations at any point in the radiation era is given by 

VA{k;v) = ^T\k;r,) (^A^ p-(fc). (343) 

Using Equation (|342p . we see that the RS matter power spectrum (evaluated in the low-energy 
regime) is ~ 50 times bigger than the GR prediction on scales given by A: ~ 10'^ /cc- 

The amplitude enhancement of perturbations is important on comoving scales < 10 AU, which 
are far too small to be relevant to present-day/cosmic microwave background measurements of 
the matter power spectrum. However, it may have an important bearing on the formation of 
compact objects such as primordial black holes and boson stars at very high energies, i.e. the 
greater gravitational force of attraction in the early universe will create more of these objects 
than in GR (different aspects of primordial black holes in RS cosmology in the context of various 
effective theories have been considered in [Ml IIMI M\ IMZl [3M1 IMS). 
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relative wavenumber log(fc/fcc) relative wavenumber log(fc/fcc) 



Figure 12: Density perturbation enhancement factors (left) and transfer functions (right) from 
simulations, effective theory, and general relativity. All of the Q factors monotonically increase 
with k/kc, and we see that the A amplitude enhancement due to 0{p/\) effects Q„t is generally 
larger than the enhancement due to KK effects Qc. For asymptotically small scales fc ^ /cc, 
the enhancement seems to level off. The transfer functions in the right panel are evaluated at a 
given subcritical epoch in the radiation dominated era. The T functions show how, for a fixed 
primordial spectrum of curvature perturbations V^^\ the effective theory predicts excess power 
in the A spectrum Va oc T'^'P'^' on supercritical/subhorizon scales compared to the GR result. 
The excess small-scale power is even greater when KK modes are taken into account, as shown by 
T^u{k',r]). Figure taken from |69| . 
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6.6 Vector perturbations 

The vorticity propagation equation on the brane is the same as in general relativity, 



UJ 



-curl Af,. 



(344) 



Taking the curl of the conservation equation (|112p (for the case of a perfect fluid, = = tt^„), 
and using the identity in Equation pSSp . one obtains 



curl A,, 



(345) 



as in general relativity, so that Equation p44|) becomes 



2-3cf Hu„ 



0, 



(346) 



which expresses the conservation of angular momentum. In general relativity, vector perturbations 
vanish when the vorticity is zero. By contrast, in brane- world cosmology, bulk KK effects can source 
vector perturbations even in the absence of vorticity [306]. This can be seen via the divergence 
equation for the magnetic part iJ^^ of the 4D Weyl tensor on the brane. 



4 2 1 2 i-£ 
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where iJ^^ = ^ (^^H^y Even when = 0, there is a source for gravimagnctic terms on the brane 
from the KK quantity curl . 

We define covariant dimensionless vector perturbation quantities for the vorticity and the KK 
gravi-vector term: 



(3f, = - curl q 



(348) 



On large scales, we can find a closed system for these vector perturbations on the brane |306) : 
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Thus we can solve for and /3p on super-Hubble scales, as for density perturbations. Vorticity in 
the brane matter is a source for the KK vector perturbation on large scales. Vorticity decays 
unless the matter is ultra-relativistic or stiffer {w > ^), and this source term typically provides a 
decaying mode. There is another pure KK mode, independent of vorticity, but this mode decays 
like vorticity. For w = p/ p = const., the solutions are 
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where 6^ = = c^. 

Inflation will redshift away the vorticity and the KK mode. Indeed, the massive KK vector 
modes are not excited during slow- roll inflation [SH] 1365] . 
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6.7 Tensor perturbations 

The covariant description of tensor modes on the brane is via the shear, which satisfies the wave 
equation |306J 



5Han„ - 



.-2 (i^£ 

1^ 



(353) 



Unlike the density and vector perturbations, there is no closed system on the brane for large scales. 
The KK anisotropic stress tt^^ is an unavoidable source for tensor modes on the brane. Thus it is 
necessary to use the 5D metric-based formalism. This is the subject of the next Section [71 
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7 Gravitational Wave Perturbations in Brane- World Cos- 
mology 

7.1 Analytical approaches 

The tensor perturbations are given by Equation (|273p . i.e., (for a flat background brane), 

(5)ds2 ^ -N'\t, y)di^ + A^it, y) [J„ + f^j] dx'dx^ + dy\ (354) 

The transverse traceless fij satisfies Equation ()293p . which implies, on splitting into Fourier 
modes with amplitude f{t,y), 
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By the transverse traceless part of Equation (|279p . the boundary condition is 

f'l 

I branc 



(355) 
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where 7fy is the tensor part of the anisotropic stress of matter-radiation on the brane. 

The wave equation (|355p cannot be solved analytically except if the background metric functions 
are separable, and this only happens for maximally symmetric branes, i.e., branes with constant 
Hubble rate Hq. This includes the RS case Hq = already treated in Section[2] The cosmologically 
relevant case is the de Sitter brane, Hq > 0. We can calculate the spectrum of gravitational 
waves generated during brane inflation [274| 11771 1149| 1233) . if we approximate slow- roll inflation 
by a succession of de Sitter phases. The metric for a de Sitter brane dS4 in AdSs is given by 
Equations (HHSl [Wl [M]) with 



7V(i,2/) = n(y), 
A{t,y)^a{t)n{y), 
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n{y) = cosh fiy ^ (^1 + y ) sinh/x|?/|, 
a{t) = aoexpiJo(i - io), 
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where fi = £~^. 

The linearized wave equation p55p is separable. As before, we separate the amplitude as 
/ = J2 fmit) fmiy) where m is the 4D mass, and this leads to: 
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The general solutions for m > are 
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Figure 13: Graviton "volcano" potential around the dS4 brane, showing the mass gap. (Figure 
taken from [272] .) 

where Bi, is a linear combination of Besscl functions, is a linear combination of associated 
Legendre functions, and 

It is more useful to reformulate Equation p63p as a Schrodinger-type equation. 



using the eonformal coordinate 



V{z)^m = -m'^m, (367) 



z = Zh+ / — — , — smh — , (368) 

Jo n{y) Ho \ J 

and defining "ifm = n'^^^fm- The potential is given by (see Figure [T5|) 

4 smh (iio^) 4 V A / 

where the last term incorporates the boundary condition at the brane. The "volcano" shape of 
the potential shows how the 5D graviton is localized at the brane at low energies. (Note that 
localization fails for an AdS4 brane |226l 1399] .) 

The non-zero value of the Hubble parameter implies the existence of a mass gap |155) , 

Am = ^Ho, (370) 

between the zero mode and the continuum of massive KK modes. This result has been generalized: 
For dS4 brane (s) with bulk scalar field, a universal lower bound on the mass gap of the KK tower 
is [149] 

Am > J^Ho. (371) 
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The massive modes decay during inflation, according to Equation p64p . leaving only the zero mode, 
which is effectively a 4D gravitational wave. The zero mode, satisfying the boundary condition 



is given by 

where the normalization condition 
implies that the function F is given by |274) 



/o(a;,0)=0, 
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At low energies {Hq <C /i) we recover the general relativity amplitude: F ^ 1. At high energies, 
the amplitude is considerably enhanced: 



Ho > At 



F ; 



2/i 



(376) 



The factor F determines the modification of the gravitational wave amplitude relative to the 
standard 4D result: 

Af2 \ 2tt 



F^{H„/^l). 



(377) 



The modifying factor F can also be interpreted as a change in the effective Planck mass |149j . 

This enhanced zero mode produced by brane inflation remains frozen outside the Hubble ra- 
dius, as in general relativity, but when it re-enters the Hubble radius during radiation or matter 
domination, it will no longer be separated from the massive modes, since H will not be constant. 
Instead, massive modes will be excited during re-entry. In other words, energy will be lost from the 
zero mode as 5D gravitons are emitted into the bulk, i.e., as massive modes are produced on the 
brane. A phenomenological model of the damping of the zero mode due to 5D graviton emission 
is given in [283] . Self-consistent low-energy approximations to compute this effect are developed 

in una ESS- 

At zero order, the low-energy approximation is based on the following [3231 13251 [27] . In the 
radiation era, at low energy, the background metric functions obey 



(378) 



To lowest order, the wave equation therefore separates, and the mode functions can be found 
analytically |323l 13251 [?7] . The massive modes in the bulk, fm{y), are the same as for a Minkowski 
brane. On large scales, or at late times, the mode functions on the brane are given in conformal 
time by 

/™„2 ,, \ 

(379) 



where ah marks the start of the low-energy regime (ph = A), and B^, denotes a linear combination 
of Bcsscl functions. The massive modes decay on super-Hubble scales, unlike the zero-mode. 
Expanding the wave equation in pa/X, one arrives at the first order, where mode-mixing arises. 
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The massive modes (pm^^ (rj) on sub- Hubble scales are sourccd by the initial zero mode that is 
re-entering the Hubble radius |137) : 

[W ~ "^"^ ^^'"^"^ ^ ^ -4P^I„.okWo \ (380) 

where Imo is a transfer matrix coefficient. The numerical integration of the equations [202] confirms 
the effect of massive mode generation and consequent damping of the zero-mode. 



7.2 Full numerical solutions 

Full numerical solutions for the tensor perturbations have been obtained by the two methods - the 
pseudo-spectral (PS) method used in |201l 1200] and the characteristic integration (CI) algorithm 
developed in [70]. It was shown that both methods give identical results and the behaviour of 
gravitational waves on a brane is quite insensitive to the initial conditions in the bulk as for the 
scalar perturbations |381j . Here we summarize the results obtained in j200) . It is convenient to use 
the static bulk metric and consider a moving brane. The simplest initial condition is f{z, t) = const 
when the mode is outside the bubble horizon on the brane. If the brane is static, this would give a 
zero-mode solution / = cos(fc(T — tq)) where tq is an initial time. However, due to the motion of the 
brane, which causes the expansion of the brane universe, KK modes are excited and this solution is 
modified. Figure [T4l demonstrate this effect. Once the perturbation enters the horizon, non-trivial 
waves are excited in the bulk and the amplitude of the tensor perturbation is damped. Figure 1151 
shows the behaviour of gravitational waves for two different wave numbers. Here = p/A at the 
time when the mode re-enters the horizon. As for scalar perturbations, we can define the effective 
4D solutions by ignoring the bulk as a reference 

h,,[ + 3Hh,,i + /ircf = 0. (381) 

a{ty 

h^ci only takes into account the effect of the high-energy modification of the Fricdmann equation. 
Figure [TS] shows that the full solution has an additional suppression of the amplitude compared 
with href. This suppression is caused by the excitations of KK modes at the horizon crossing as is 
seen in Figure 1141 The suppression is stronger for modes that enter the horizon earlier e, > 1 and 
it becomes negligible at low energies e, < 1. 

The ratio |/i5D//ircf| evaluated at the low-energy regime long after the horizon re-entry time 
monotonically decreases with the frequency and the suppression of amplitude /15D becomes signif- 
icant above the critical frequency /crit given by 
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(382) 



This corresponds to a frequency of the mode that enters the horizon when iJ* = £ i (cf. |203) ). 
The ratio |/i5D/^rof| obtained from numerical solutions is fitted as 



«5D 



/crit 



(383) 



with a = 0.76 ± 0.01 and /3 = 0.67 ± 0.01. 
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Figure 14: The evolution of gravitational waves. We set the comoving wave number to fc = ^/2l/^ 
(e* = 1.0). The right panel depicts the projection of the three-dimensional waves of the left panel. 
Figure taken from |200| . 




Figure 15: Squared amplitude of gravitational waves on the brane in the low-energy {left) and 
the high-energy {right) regimes. In both panels, solid lines represent the numerical solutions. The 
dashed lines are the amplitudes of reference gravitational waves /ircf obtained from Equation p81|) . 
Figure taken from |200j . 
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There are two important effects on the spectrum in the high-energy regime. Let us first consider 
the non-standard cosmological expansion due to the p^-tcrm. The spectrum of the stochastic 
gravitational waves is modified to 

^..4(2 f (384) 

where w is an equation of state. This is because gravitational waves re-enter the horizon when the 
p^-term dominates at high frequencies / > /crit- In the high-energy radiation dominated phase, 
the spectrum of the stochastic gravitational waves is modified to 

f^rcf « (/crit < /). (385) 

The other effect is the KK-modc excitations. Taking account of the KK-mode excitations, the 
spectrum is calculated as 



GW 



Kd 



^rof 



nrel, (386) 



where we used the fact JIqw oc h? P . Combining it with the result (|383p . the spectrum becomes 
nearly flat above the critical frequency: 

^GW « /°, (387) 

which is shown in filled squares in Figure 1161 In this figure, the spectrum calculated from the 
reference gravitational waves ilrof is also shown in filled circles. Note that the normalization factor 
of the spectrum is determined as JIqw = 10~^^ from the CMB constraint. The short-dashed line 
and the solid line represent asymptotic behaviors in the high-frequency region. The spectrum 
taking account of the two high-energy effects seems almost indistinguishable from the standard 
four-dimensional prediction shown in long-dashed line in the figure. In other words, while the 
effect due to the non-standard cosmological expansion enhances the spectrum, the KK-mode effect 
reduces the GW amplitude, which results in the same spectrum as the one predicted in the four- 
dimensional theory. Note that the amplitude taking account of the two effects near / w /crit is 
slightly suppressed, which agrees with the results in the previous study for < 0.3 using the 
Gaussian-normal coordinates |202) discussed in the previous subsection. 

This cancelation of two high energy effects is valid only for w = 1/3. For other equations of 
state, the final spectrum at high frequencies are different from 4D predictions. For example for 
u) = 1, f^GW p/^ for f > fc while the 4D theory predicts flow oc p. 
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Figure 16: The energy spectrum of the stochastic background of gravitational wave around the 
critical frequency fc in radiation dominated epoch. The filled circles represent the spectrum caused 
by the non-standard cosmological expansion of the universe. Taking account of the KK-mode 
excitations, the spectrum becomes the one plotted by filled squares. In the asymptotic region 
depicted in the solid line, the frequency dependence becomes almost the same as the one predicted 
in the four-dimensional theory (long-dashed line). Figure taken from [200j . 

8 CMB Anisotropies in Brane- World Cosmology 

For the CMB anisotropies, one needs to consider a multi-component source. Linearizing the general 
nonlinear expressions for the total effective energy-momentum tensor, we obtain 



Ptot =P+^i2p + p) + y, 



1 - 



p + 'ip 

2A 



where 



(388) 
(389) 
(390) 

(391) 
(392) 



are the total matter-radiation density, pressure, and momentum density, respectively, and tt^i, is 
the photon anisotropic stress (neglecting that of neutrinos, baryons, and CDM). 

The perturbation equations in the previous Section [7] form the basis for an analysis of scalar and 
tensor CMB anisotropies in the brane-world. The full system of equations on the brane, including 
the Boltzmann equation for photons, has been given for scalar j284| and tensor [283] perturbations. 
But the systems are not closed, as discussed above, because of the presence of the KK anisotropic 
stress TT^j^, which acts a source term. 
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In the tight-couphng radiation era, the scalar perturbation equations may be decoupled to give 
an equation for the gravitational potential defined by the electric part of the brane Wcyl tensor 
(not to be confused with £^^): 

= V(^V,)$. (393) 

In general relativity, the equation in $ has no source term, but in the brane- world there is a source 
term made up of tt^^ and its time-derivatives. At low energies (p <C A), and for a flat background 
(A' = 0), the equation is [284] 



const. 
A 



(394) 



where x ~ k/{aH), a prime denotes d/dx, and and tt^ are the Fourier modes of $ and tt^^, 
respectively. In general relativity the right hand side is zero, so that the equation may be solved 
for and then for the remaining perturbative variables, which gives the basis for initializing 
CMB numerical integrations. At high energies, earlier in the radiation era, the decoupled equation 
is fourth order |284| : 

729x'^<i>'l" + 3888x^';[' + (1782 + Mx^)<i>k + 1442;$'^ + (90 + x^)<i>k 
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(395) 



The formalism and machinery are ready to compute the temperature and polarization anisotro- 
pics in branc-world cosmology, once a solution, or at least an approximation, is given for tt^^. The 
resulting power spectra will reveal the nature of the brane- world imprint on CMB anisotropics, and 
would in principle provide a means of constraining or possibly falsifying the branc-world models. 
Once this is achieved, the implications for the fundamental underlying theory, i.e., M theory, would 
need to be explored. 

However, the first step required is the solution for tt^^,. This solution will be of the form 
given in Equation (|249p . Once Q and Fk are determined or estimated, the numerical integration 
in Equation p49p can in principle be incorporated into a modified version of a CMB numerical 
code. The full solution in this form represents a formidable problem, and one is led to look for 
approximations . 



8.1 The low-energy approximation 



The basic idea of the low-energy approximation j401l 14311 13901 14021 I403J is to use a gradient 
expansion to exploit the fact that, during most of the history of the universe, the curvature scale 
on the observable brane is much greater than the curvature scale of the bulk {£ < 1 mm) : 



These conditions are equivalent to the low energy regime, since 
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Using Equation p96p to neglect appropriate gradient terms in an expansion in £^/L^, the low- 
energy equations can be solved. However, two boundary conditions are needed to determine all 
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functions of integration. This is achieved by introducing a second brane, as in the RS 2-brane 
scenario. This brane is to be thought of cither as a regulator brane, whose backreaction on the 
observable brane is neglected (which will only be true for a limited time), or as a shadow brane 
with physical fields, which have a gravitational effect on the observable brane. 

The background is given by low-energy FRW branes with tensions ±A, proper times t±, scale 
factors a± , energy densities p± and pressures p± , and dark radiation densities p£ ± . The physical 
distance between the brancs is £d{t), and 
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Then the background dynamics is given by 
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(399) 
(400) 



(see [39j 1278) for the general background, including the high-energy regime). The dark energy 
obeys p5+ = C/a^, where C is a constant. From now on, we drop the -I— subscripts which refer to 
the physical, observed quantities. 

The perturbed metric on the observable (positive tension) brane is described, in longitudinal 
gauge, by the metric perturbations ip and TZ, and the perturbed radion is c? = d + N. The 
approximation for the KK (Weyl) energy-momentum tensor on the observable brane is 
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and the field equations on the observable brane can be written in scalar-tensor form as 
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where 
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The perturbation equations can then be derived as generalizations of the standard equations. 
For example, the SG'^q equation is |245| 



1 k'^ 

H^ij -HTZ 

3 a-' 



1 



-K 



.2 



^2d 



6 e 
1 



2d 



5p + e '^'^Sp- ) + -K 



^2d 



a2d 



d- H] N 



3 e 

2 



2d _ I 
J2 



PsN 



d-H] N - d^il! + 2Hd%i^- dn- -—N 



IP 
3a2 



(404) 



The trace part of the perturbed field equation shows that the radion perturbation determines the 
crucial quantity, 5tt£: 



:2d , 



1 



(405) 
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where the last equahty foUows from Equation p25p . The radion perturbation itself satisfies the 
wave equation 



N + [3H-2d 



N- (2H + + 2d-' - 6Hd - 2d] N 
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A new set of variables (p± , E turns out be very useful |246l 1247] : 
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Equation (|405p gives 
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The variable E determines the metric shear in the bulk, whereas ip± give the brane displace- 
ments in transverse traceless gauge. The latter variables have a simple relation to the curvature 
perturbations on large scales |246l I247J (restoring the +- subscripts) : 



(409) 



where f± — df±/dt±. 



8.2 The simplest model 

The simplest model is the one in which 

P£=0 = d (410) 

in the background, with p^/p^ = p/p. The regulator brane is assumed to be far enough away that 
its effects on the physical brane can be neglected over the timescales of interest. By Equation (|400p 
it follows that 

P- = -pe'', (411) 

i.e., the matter on the regulator brane must have fine-tuned and negative energy density to prevent 
the regulator brane from moving in the background. With these assumptions, and further assuming 
adiabatic perturbations for the matter, there is only one independent brane- world parameter, i.e., 
the parameter measuring dark radiation fluctuations: 

SC, ^ (412) 

Prad 

This assumption has a remarkable consequence on large scales: The Weyl anisotropic stress 
Sirs terms in the Sachs-Wolfe formula p27p cancel the entropy perturbation from dark radiation 
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Figure 17: The CMB power spectrum with brane-world effects, encoded in the dark radiation 
fluctuation parameter SC^ as a proportion of the large-scale curvature perturbation for matter 
(denoted C* hi the plot). (Figure taken from [247] .) 



fluctuations, so that there is no difference on the largest scales from the standard general relativity 
power spectrum. On small scales, beyond the first acoustic peak, the branc-world corrections arc 
negligible. On scales up to the first acoustic peak, brane-world effects can be significant, changing 
the height and the location of the first peak. These features are apparent in Figure[T71 However, it 
is not clear to what extent these features are general brane-world features (within the low-energy 
approximation), and to what extent they are consequences of the simple assumptions imposed on 
the background. Further work remains to be done. 

A related low-energy approximation, using the moduli space approximation, has been developed 
for certain 2-branc models with bulk scalar field |364[ \5T\ . The effective gravitational action on 
the physical branc, in the Einstein frame, is 



5cff = ^-^ / d'^x 
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where a is a coupling constant, and 
the bulk scalar field and the radion). 
the x-field. 



(j) and X s-rc moduli fields (determined by the zero-mode of 
Figure [18] shows how the CMB anisotropics are affected by 
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9 DGP Models: Modifying Gravity at Low Energies 



9.1 'Self-accelerating' DGP 

Most branc- world models modify general relativity at high energies. The Randall-Sundrum models 
discussed up to now are a typical example. At low energies, Hi <C 1, the zero- mode of the graviton 
dominates on the brane, and general relativity is recovered to a good approximation. At high 
energies, H(. 3> 1, the massive modes of the graviton dominate over the zero-mode, and gravity 
on the brane behaves increasingly 5-dimensional. On the unperturbed FRW brane, the standard 
energy-conservation equation holds, but the Friedmann equation is modified by an ultraviolet 
correction, {Gipf'. At high energies, gravity "leaks" off the brane and H'^ oc . 

By contrast, the brane-world model of Dvali-Gabadadze-Porrati |135| (DGP), which was gen- 
eralized to cosmology by Deffayet |116) . modifies general relativity at low energies. This model 
produces 'self-acceleration' of the late-time universe due to a weakening of gravity at low energies. 
Like the RS model, the DGP model is a 5D model with infinite extra dimensions!! 

The action is given by 



167rG 



J bulk J branc 



(414) 



The bulk is assumed to be 5D Minkowski spacetime. Unlike the AdS bulk of the RS model, the 
Minkowski bulk has infinite volume. Consequently, there is no normalizable zero-mode of the 4D 
graviton in the DGP brane-world. Gravity leaks off the 4D brane into the bulk at large scales, 
r ^ Tc, where the first term in the sum (|414p dominates. On small scales, gravity is effectively 
bound to the brane and 4D dynamics is recovered to a good approximation, as the second term 
dominates. The transition from 4D to 5D behaviour is governed by the crossover scale re- For a 
Minkowski brane, the weak-field gravitational potential behaves as 

r,-ifo,,«,^ . (415) 

yr ^ toT r :$> Vc 

On a Friedmann brane, gravity leakage at late times in the cosmological evolution can initiate 
acceleration - not due to any negative pressure field, but due to the weakening of gravity on the 
brane. 4D gravity is recovered at high energy via the lightest massive modes of the 5D graviton, 
effectively via an ultra-light meta-stable graviton. 

The energy conservation equation remains the same as in general relativity, but the Friedmann 
equation is modified I116j : 

p + 3H{p + p) = 0, (416) 

The Friedmann equation can be derived from the junction condition or the Gauss-Codazzi equation 
as in the RS model. To arrive at Equation (|417l) we have to take a square root which implies a 
choice of sign. As we shall see, the above choice has the advantage of leading to acceleration - 
but at the expense of a 'ghost' (negative kinetic energy) mode in the scalar graviton sector. The 
'normal' (non-self-accelerating) DGP model, where the opposite sign of the square root is chosen, 
has no ghost, and is discussed below. 

From Equation ()417|) we infer that at early times, i.e., Hrc ^ 1, the general relativistic 
Friedman equation is recovered. By contrast, at late times in an expanding CDM universe, with 



^DGP brane-worlds without Z2 symmetry have also been considered; see e.g. |396| . 
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Figure 19: Joint constraints [solid thick (blue)] from the SNLS data [solid thin (yellow)], the 
BAO peak at z = 0.35 [dotted (green)] and the CMB shift parameter from WMAP3 [dot-dashed 
(red)]. The left plot shows LCDM models, the right plot shows DGP. The thick dashed (black) 
line represents the flat models, Hk — 0. (From [308].) 




Figure 20: The constraints from SNe and CMB/BAO on the parameters in the DGP model. The 
flat DGP model is indicated by the vertical dashed-dotted line; for the MFCS light curve fit, the 
flat model matches to the date very well. The SAFT-II light curve fit to the SNe is again shown 
by the dotted contours. The combined constraints using the SAFT-II SNe outlined by the dashed 
contours represent a poorer match to the CMB/BAO for the flat model. (From |404| .) 
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p (X a — > 0, we have 

H^Hoo^-. (418) 
so that expansion accelerates and is asymptoticaUy de Sitter. The above equations imply 
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H- — ^ -AnGp 



1 



(419) 



^1 + 327rGr2p/3 
In order to achieve self-acceleration at late times, we require 

> H^' , (420) 



since Hq ^ Hoc- This is confirmed by fitting supernova observations, as shown in Figures [19] and 
1201 The dimensionless cross-over parameter is defined as 

and the LCDM relation, 

+ Oa + t^A' = 1 , (422) 

is modified to 



Qm, + 2ya~Vl - ^^A' + Oa' = 1 . (423) 

LCDM and DGP can both account for the supernova observations, with the fine-tuned values 
A ^ Hq and rc ^ Hq^, respectively. When we add further constraints to the expansion history 
from the baryon acoustic oscillation peak at z = 0.35 and the CMB shift parameter, the DGP fiat 
models are in strong tension with the data, whereas LCDM models provide a consistent fit. This 
is evident in Figures [19] and [20] though this conclusion depends on a choice of light curve fitters 
in the SNe observations. The open DGP models provide a somewhat better fit to the geometric 
data - essentially because the lower value of r2,„ favoured by supernovae reduces the distance to 
last scattering and an open geometry is able to extend that distance. For a combination of SNe, 
CMB shift and Hubble Key Project data, the best-fit open DGP also performs better than the fiat 
DGP gog, as shown in Figure [2l] 

Observations based on structure formation provide further evidence of the difference between 
DGP and LCDM, since the two models suppress the growth of density perturbations in different 
ways |2981 1297] . The distance-based observations draw only upon the background 4D Friedman 
equation (|417p in DGP models - and therefore there are quintessence models in general relativity 
that can produce precisely the same supernova distances as DGP. By contrast, structure formation 
observations require the 5D perturbations in DGP, and one cannot find equivalent quintessence 
models |253) . One can find 4D general relativity models, with dark energy that has anisotropic 
stress and variable sound speed, which can in principle mimic DGP |265) . However, these models 
are highly unphysical and can probably be discounted on grounds of theoretical consistency. 

For LCDM, the analysis of density perturbations is well understood. For DGP the perturbations 
are much more subtle and complicated j253j . Although matter is confined to the 4D brane, gravity 
is fundamentally 5D, and the 5D bulk gravitational field responds to and back-reacts on 4D density 
perturbations. The evolution of density perturbations requires an analysis based on the 5D nature 
of gravity. In particular, the 5D gravitational field produces an effective "dark" anisotropic stress 
on the 4D universe, as discussed in Section 13.41 If one neglects this stress and other 5D effects, 
and simply treats the perturbations as 4D perturbations with a modified background Hubble rate 
- then as a consequence, the 4D Bianchi identity on the brane is violated, i.e., V^G^^ ^ 0, and the 
results are inconsistent. When the 5D effects are incorporated |253l 171) . the 4D Bianchi identity 
is automatically satisfied. (See Figure]^) 

There are three regimes governing structure formation in DGP models: 
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Figure 21: The difference in between best- fit DGP (flat and open) and best-fit (flat) LCDM, 
using SNc, CMB shift and Hq Key Project data. (From gDS].) 
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• On small scales, below the Vainshtein radius (which for cosmological purposes is roughly the 
scale of clusters), the spin-0 scalar degree of freedom becomes strongly coupled, so that the 
general relativistic limit is recovered [258) . 

• On scales relevant for structure formation, i.e., between cluster scales and the Hubble ra- 
dius, the spin-0 scalar degree of freedom produces a scalar-tensor behaviour. A quasi-static 
approximation (as in the Newtonian approximation in standard 4D cosmology) to the 5D 
perturbations shows that DGP gravity is like a Brans-Dicke theory with parameter |253j 

UJBD = ^{I3-1), 

13 = 1 + 2H^r, (^H^ + J 

At late times in an expanding universe, when Hrc ^ 1, it follows that /3 < 1, so that ujbd < 0. 
(This is a signal of the ghost pathology in DGP, which is discussed below.) 

• Although the quasi-static approximation allows us to analytically solve the 5D wave equation 
for the bulk degree of freedom, this approximation breaks down near and beyond the Hubble 
radius. On super-horizon scales, 5D gravity effects arc dominant, and we need to solve 
numerically the partial differential equation governing the 5D bulk variable [7T] . 
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Figure 22: The growth factor g{a) = A{a)/a for LCDM (long dashed) and DGP (sohd, thick), as 
well as for a dark energy model with the same expansion history as DGP (short dashed). DGP-4D 
(solid, thin) shows the incorrect result in which the 5D effects are set to zero. (From |253j .) 

On subhorizon scales relevant for linear structure formation, 5D effects produce a difference 
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Figure 23: Left: Numerical solutions for DGP density and metric perturbations, showing also 
the quasi-static solution, which is an increasingly poor approximation as the scale is increased. 
(From [7T].) Right: Constraints on DGP (the open model in Figure [5T] that provides a best fit to 
geometric data) from CMB anisotropies (WMAP5). The DGP model is the sohd curve, QCDM 
(short-dashed curve) is the GR quintessence model with the same background expansion history 
as the DGP model, and LCDM is the dashed curve (a slightly closed model that gives the best fit 
to WMAP5, HST and SNLS data). (From [143].) 
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between (p and — |253j : 



3/3 

so that there is an effective dark anisotropic stress on the brane: 



AnGa^l-j^^pA, (426) 
/c^V = -47rGa^ ( 1 + ) pA , (427) 



fc2(0 + ^) = -^^pA. (428) 



The density perturbations evolve as 

A + 2i/A -47rG ^1 - pA = 0. (429) 

The linear growth factor, g{a) = A(a)/a (i.e., normalized to the flat CDM case, A cx a), is shown 
in Figure [22] This illustrates the dramatic suppression of growth in DGP relative to LCDM - 
from both the background expansion and the metric perturbations. If we parametrize the growth 
factor in the usual way, we can quantify the deviation from general relativity with smooth dark 
energy [295] : 

dlnA ^ , ,^ f 0.55 + 0.05(1 + u>(z = 1)1 GR, smooth DE 

dh^ = "'"('^^ ' 10.68 DGP (4^0) 

Observational data on the growth factor [189] are not yet precise enough to provide meaningful 
constraints on the DGP model. Instead, we can look at the large-angle anisotropics of the CMB, 
i.e., the ISW effect. This requires a treatment of perturbations near and beyond the horizon 
scale. The full numerical solution has been given by [7T], and is illustrated in Figure [23] The 
CMB anisotropics are also shown in Figure [23] as computed in |143j using a scaling approximation 
to the supcr-Hubble modes j379) (the accuracy of the scaling ansatz is verified by the numerical 
results f7T\). 

It is evident from Figure that the DGP model, which provides a best fit to the geometric 
data (see Figure I^Tj) . is in serious tension with the WMAP5 data on large scales. The problem 
arises because there is a large deviation of 0- = (0 — ijj)/2 in the DGP model from the LCDM 
model. This deviation, i.e., a stronger decay of (/)_, leads to an over-strong ISW effect (which is 
determined by t/)_), in tension with WMAP5 observations. 

As a result of the combined observations of background expansion history and large-angle 
CMB anisotropics, the DGP model provides a worse fit to the data than LCDM at about the 5cr 
level |143j . Effectively, the DGP model is ruled out by observations in comparison with the LCDM 
model. 

In addition to the severe problems posed by cosmological observations, a problem of theo- 
retical consistency arises from the fact that the late-time asymptotic dc Sitter solution in DGP 
cosmological models has a ghost. The ghost is signaled by the negative Brans-Dicke parameter 
in the effective theory that approximates the DGP on cosmological subhorizon scales: The ex- 
istence of the ghost is confirmed by detailed analysis of the 5D perturbations in the dc Sitter 
limit [2481 11761 ISTl 1248] . The DGP ghost is a ghost mode in the scalar sector of the gravitational 
field - which is more serious than the ghost in a phantom scalar field. It effectively rules out the 
DGP, since it is hard to see how an ultraviolet completion of the DGP can cure the infrared ghost 
problem. However, the DGP remains a valuable toy model for illustrating the kinds of behaviour 
that can occur from a modification to Einstein's equations - and for developing cosmological tools 
to test modified gravity and Einstein's theory itself. 
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9.2 'Normal' DGP 



The self-accelerating DGP is effectively ruled out as a viable cosmological model by observations 
and by the problem of the ghost in the gravitational sector. Indeed, it may be the case that 
self-acceleration gencrically comes with the price of ghost states. The 'normal' (i.e., non-sclf- 
accelerating and ghost- free) branch of the DGP j374) . arises from a different embedding of the 
DGP brane in the Minkowski bulk (see Figure [M]). In the background dynamics, this amounts 
to a replacement Vc — > —Vc in Equation (|417p - and there is no longer late-time self-acceleration. 
Therefore, it is necessary to include a A term in order to accelerate the late universe: 

Tr2 K 1 ri F SttG a 

(Normal DGP models with a quintessence field have also been investigated [HH]-) Using the di- 
mensionlcss crossover parameter defined in Equation (|42ip . the densities are related at the present 
time by 

v/l - fix = - + V^^r, + a„ + f^A , (432) 

which can be compared with the self-accelerating DGP relation (j423|) . 




Figure 24: Left: The embedding of the self-accelerating and normal branches of the DGP brane 
in a Minkowski bulk. (From [81].) Right: Joint constraints on normal DGP (fiat, K = Q) from 
SNLS, CMB shift (WMAP3) and BAG {z = 0.35) data. The best-fit is the solid point, and is 
indistinguishable from the LGDM limit. The shaded region is unphysical and its upper boundary 
represents fiat LCDM models. (From [281] .) 

An interesting feature of the normal branch is the 'degravitation' property, i.e., that A is 
effectively screened by 5D gravity effects. This follows from rewriting the modified Friedmann 
equation (|43ip in standard general relativistic form, with 

Ks = K^-\ m + \<k. (433) 
Tc V a 

Thus, 5D gravity in normal DGP can in principle reduce the bare vacuum energy significantly |374[ 
1300) . However, Figure shows that best-fit fiat models, using geometric data, only admit insignif- 
icant screening [281) . The closed models provide a better fit to the data [168) . and can allow a 
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bare vacuum energy term with fl\ > 1, as shown in Figure [25l This does not address the funda- 
mental problem of the smallness of 51a, but is nevertheless an interesting feature. We can define 
an effective equation-of-state parameter via 

Aeff + 3H{1 + Wcff )Aeff = . (434) 

At the present time (setting K — for simplicity), 

"^•^^^ = - (i-a.)(a„ + fiA + i) < ' ^'''^ 

where the inequality holds because 57^ < 1- This reveals another important property of the nor- 
mal DGP model: effective phantom behaviour of the recent expansion history. This is achieved 
without any pathological phantom field (similar to what can be done in scalar-tensor theories [42]). 
Furthermore, there is no "big rip" singularity in the future associated with this phantom acceler- 
ation, unlike the situation that typically arises with phantom fields. The phantom behaviour in 
the normal DGP model is also not associated with any ghost problem ~ indeed, the normal DGP 
branch is free of the ghost that plagues the self-accelerating DGP [HI] ■ 
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Figure 25: Left: Joint constraints on normal DGP from SNc Gold, CMB shift (WMAP3) and Hq 
data in the projected curvature-A plane, after marginalizing over other parameters. The best-fits 
are the solid points, corresponding to different values of ilm- (From |168) .) Right: Numerical 
solutions for the normal DGP density and metric perturbations, showing also the quasistatic so- 
lution, which is an increasingly poor approximation as the scale is increased. Compare with the 
self-accelerating DGP case in Figure (From [71].) 
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Perturbations in the normal branch have the same structure as those in the self-accelerating 
branch, with the same regimes - i.e., below the Vainshtein radius (recovering a GR limit), up to the 
Hubble radius (Brans-Dickc behaviour), and beyond the Hubble radius (strongly 5D behaviour). 
The quasistatic approximation and the numerical integrations can be simply repeated with the 
replacement Vc —Tc (and the addition of A to the background). In the sub- Hubble regime, 
the effective Brans-Dicke parameter is still given by Equations (|424|) and (|425|) . but now we have 
uiBD > - and this is consistent with the absence of a ghost. Furthermore, a positive Brans-Dicke 
parameter signals an extra positive contribution to structure formation from the scalar degree of 
freedom, so that there is less suppression of structure formation than in LCDM - the reverse of 
what happens in the self-accelerating DGP. This is confirmed by computations, as illustrated in 
Figure [211 

The closed normal DGP models fit the background expansion data reasonably well, as shown 
in Figure 1251 The key remaining question is how well do these models fit the large-angle CMB 
anisotropics, which is yet to be computed at the time of writing. The derivative of the ISW potential 
0_ can be seen in Figure [25l and it is evident that the ISW contribution is negative relative to 
LCDM at high redshifts, and goes through zero at some rcdshift before becoming positive. This 
distinctive behaviour may be contrasted with the behaviour in f{R) models (see Figure [26]) : both 
types of model lead to less suppression of structure than LCDM, but they produce different ISW 
effects. However, in the limit — >■ c«, normal DGP tends to ordinary LCDM, hence observations 
which fit LCDM will always just provide a lower limit for re- 
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Figure 26: Top: Measurement of the cross-correlation functions between six different galaxy data 
sets and the CMB, reproduced from [167] . The curves show the theoretical predictions for the 
ISW-galaxy correlations at each rcdshift for the LCDM model (black, dashed) and the three nDGP 
models, which describe the 1 — region of the geometry test from Figure E5] (From |168| .) Bottom: 
Theoretical predictions for a family of f{R) thccggcs compared with the ISW data |167j measuring 
the angular CCF between the CMB and six galaxy catalogues. The model with i?o = is equivalent 
to LCDM, while increasing departures from GR produce negative cross-correlations. (From [166].) 



10 6-Dimensional Models 



For branc-world models in 6-dinicnsional spacctimc, the codimension of a brane is two and the 
behaviour of gravity is quahtatively very different from the codimension one brane-world models. 
Here we briefly discuss some important examples and features of 6-dimensional models. 



10.1 Supersymmetric Large Extra Dimensions (SLED) Model 

This is a supersymmetric version of the Einstein-Maxwell model [72] (see [57l [58] [59] 1249) for 
reviews). The bosonic part of the action is given by [3] [4] [2] 



S = / <fx^/^ 



.2 
^6 



(436) 



There exists a solution where the dilaton (jj is constant, (j) — (j)o. The solution for the 6D spacetime 
is given by 

ds^ ^ ij^^dx^'dx" + jijdx'^dx^. (437) 

The gauge field is taken to consist of magnetic fiux threading the extra dimensional space so that 
the field strength takes the form 

F.j^^Boe,,, (438) 

where Bq is a constant, 7 is the determinant of and is the antisymmetric tensor normalized 
as £12 = 1. All other components of Fat vanish. A static and stable solution is obtained by choosing 
the extra-dimensional space to be a two-sphere 

j.jdx'dx^ = alide^ + sin^ 9dip^). (439) 

The magnetic field strength Bq and the radius qq are fixed by the cosmological constant 

Sg = 2Aee^*o, = (440) 

The constant value 0o is determined by the condition that the potential for has minimum [154j 

V'{cj>o) = -\Ble~^« + Aee*" = 0. (441) 

This is exactly the condition to have a flat geometry on the brane (see Equation (|440p ) 

S^e~*° = 2A6e'^°. (442) 

Thus without tuning the 6D cosmological constant, it is possible to obtain a flat 4D spacetime. 
Now we add branes to this solution [72] . The brane action is given by 

j4 



54 = -cr J d^Xy^. (443) 

The solution for the extra dimensions is now given by 

-f.jdx'dx^ = alide^ + sin^ Odtp"^), {AAA) 

where 

a = l-^^, al = ^^. (445) 
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The coordinate ranges from to 2tt. Thus the effect of the brane makes a deficit angle S = 
27r(f — a) in the bulk (see Figure [57]). This is a 6D realization of the ADD model including the 
self-gravity of brancs. An interesting property of this model is that regardless of the tension of 
the brane, the 4D specetime on the branc is flat. Thus this could solve the cosmological constant 
problem - any vacuum energy on the branc only changes a geometry of the extra-dimensions and 
does not curve the 4D spacetime. This idea of solving the cosmological constant problem is known 
as self-tuning. 



Figure 27: Removing a wedge from a sphere and identifying opposite sides to obtain a rugby ball 
geometry. Two equal-tension branes with conical deficit angles are located at either pole; outside 
the branes there is constant spherical curvature. From [72] . 

Wc should note that there have been several objections to the idea of self-tuning |154| 14221 1423) . 
Consider that a phase transition occurs and the tension of the brane changes from ai to a2- 
Accordingly, a changes from ai = 1 — ai/ {2t:Mq) to a2 = 1 — cr2/(27rM|). The magnetic flux is 
conserved as the gauge field strength is a closed form, dF = 0. Then the magnetic flux which is 
obtained by integrating the field strength over the extra dimensions should be conserved 

= 47rQ:ii3o.i = 47rQ;2-Bo,2- (446) 

The relation between Ag and Bq , Ec^uation ()440p , that ensures the existence of Minkowski branes 
cannot be imposed both for Bq = Bq^i and Bq = i3o,2 when ai ^ unless </)o changes. Moreover, 
the quantization condition must be imposed on the flux $5. What happens is that a modulus, 
which is a combination of and the radion describing the size of extra-dimension, acquires a 
runaway potential and the 4D spacetime becomes non-static. 

An unambiguous way to investigate this problem is to study the dynamical solutions directly in 
the 6D spacetime. However, once we consider the case where the tension becomes time dependent, 
we encounter a difficulty to deal with the branes |422j . This is because for co-dimension 2 branes, 
we encounter a divergence of metric near the brane if we put matter other than tension on a brane. 
Hence, without specifying how we regularize the branes, we cannot address the question what 
will happen if we change the tension. Is the self-tuning mechanism at work and does it lead to 
another static solution? Or do we get a dynamical solution driven by the runaway behaviour of 
the modulus field? 

There was a negative conclusion on the self-tuning in this supersymmetric model for a particular 
kind of regularization |422] 1423] . However, the answer could depend on the regularization of 
branes and the jury remains out. It is important to study the time-dependent dynamics in the 6D 
spacetime and the regularization of the branes in detail l60 ] l62 ] lilt ] [6T ] [33 ] fTTT ] IMT ] l353 ] 
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10.2 Cosmology in 6D brane-world models 

It is much harder to obtain cosmological solutions in 6D models than in 5-dimensional models. In 
5D brane-world models, cosmological solutions can be obtained by considering a moving brane in a 
static bulk spacctime. This is because the motion of the homogeneous and isotropic brane does not 
change the bulk spacetime thanks to Birkoff 's theorem. However, this is no longer the case in 6D 
spacetime. Thus we need to find time dependent bulk solutions that are coupled to a motion of a 
brane, which requires us to solve non-linear partial differential equations numerically. Moreover, as 
is mentioned above, there appears a curvature singularity if one considers an infinitely thin brane 
and puts matter other than tension on the brane |422| 1423) . Thus we need some regularization 
scheme to find cosmological solutions. 

One of the popular ways to regularize a brane is to promote a brane that is a point-like object 
in two extra-dimensions to a 5D ring |356| . The ring is a codimension one object and it is possible 
to consider a motion of this ring. However, a problem is that this ring brane is not homogeneous 
as one of dimensions on the brane is compact and this breaks Birkhoff's theorem. It has been 
shown that cosmology obtained by the motion of a brane in a static bulk shows pathological 
behaviour [3541 1328] . This indicates that we need to find fully time-dependent bulk solutions. 

At the moment, the only accessible way is to solve the 6D bulk spacetime using the gradient 
expansion methods assuming physical scales on a brane are much lower than mass scales in the 
bulk [1511 [TBI 1234) . It is still an open question what are high energy effects in 6D models. 

As in 5D models, there are many generalizations such as the inclusion of the induced gravity 
term on a branc and the Gauss-Bonnet term in the bulk [43| [2231 HHIl EM [352] [240] [83] [82]. 



10.3 Cascading brane-world model 

This is a 6D extension of the DGP model [115[ 1112] . The action is given by 



J ^'*5 J4-bianc 



1 

'2ni 



d^x^/^i^ ^*^R + L„ 



3 — branc 



(447) 



As in the 5D DGP model, we can define the cross over scale, 
over scales: 



(5), 



(6), 



Ml 
Mi- 



In this model, there are two cross 



(448) 



Assuming that 



(6), 



> 



(5), 



the gravitational potential on the 3-brane cascades from a 1/r (4D 



gravity) regime at short scales, to a 1/r^ (5D gravity) regime at intermediate distance and finally 
a 1/r'^ (6D gravity) regime at large distances. 

This model addresses several fundamental issues in induced gravity models in 6D spacetime. 
Without the induced gravity term on the 4-brane, the 6D graviton propagator diverges logarith- 
mically near the position of the 3-brane |163) . On the other hand, the graviton propagator on 
the 3-brane in this model behaves like G{p) — >■ log(p'^^Vc) in the M5 — limit where p is 4D 
momentum. Thus the cross-over scale '•^Vc acts as a cut-off for the propagator and it remains 
finite even at the position of the 3-brane. 

A more serious issue in the induced gravity model is that most constructions seem to be plagued 
by ghost instabilities |129[I153] . Usually, the regularization of the brane is necessary and it depends 
on the regularization scheme whether there appears a ghost or not [244| . In the cascading model, 
there is still a ghost if the 3-brane has no tension. However, it was shown that by adding a tension 
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to the 3-brane, this ghost disappears |1151 11121 1113) . More precisely, there is a critical tension 



Acdt = ^^'V-2A.f| (449) 

and above the critical tension A > Adit, the model is free of ghosts. It is still unclear what is the 
meaning of having the critical tension for the existence of the ghost. For example, what happens 
if there is a phase transition on the 3-brane and the tension changes from A > Acrit to A < Acrit? 
This remains a very interesting open question. This model also provides an interesting insight 
into the "degravitation" mechanism by which the cosmological constant does not gravitate on the 
3-brane [TMlfTT^ . 

Cosmological solutions in the cascading brane model are again notoriously difficult to find be- 
cause it is necessary to find 6-dimensional solutions that depend on time and two extra-coordinates [1] . 
The simplest de Sitter solutions have been obtained |326| . Interestingly, there exists a self- 
accelerating solution for the 3-brane even when the solution for the 4-brane is in the normal 
branch. It is still not clear whether this self-accelerating solution is stable or not and it is crucial 
to check the stability of this new self-accelerating solution. 

A similar class of models includes intersecting branes |2181 11021 1103) . In this model, we have 
two 4-branes that intersect and a 3-brane sits at the intersection. Again there are self-accelerating 
de Sitter solutions and cosmology has been studied by considering a motion of one of the 4-branes. 
A model without a 4-brane has been studied by regularizing a 3-brane by promoting it to a 5D 
ring brane [HOI HE]- 
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11 Conclusion 



Simple brane-world models provide a rich phenomenology for exploring some of the ideas that 
are emerging from M theory. The higher-dimensional degrees of freedom for the gravitational 
field, and the confinement of standard model fields to the visible brane, lead to a complex but 
fascinating interplay between gravity, particle physics, and geometry, that enlarges and enriches 
general relativity in the direction of a quantum gravity theory. 

This review has attempted to show some of the key features of brane-world gravity from the 
perspective of astrophysics and cosmology, emphasizing a geometric approach to dynamics and 
perturbations. It has focused mainly on 1-brane RS-type brane-worlds, but also considered the 
DGP brane-world models. The RS-type models have some attractive features: 

• They provide a simple 5D phenomenological realization of the Horava-Witten supergravity 
solutions in the limit where the hidden brane is removed to infinity, and the moduli effects 
from the 6 further compact extra dimensions may be neglected. 

• They develop a new geometrical form of dimensional reduction based on a strongly curved 
(rather than flat) extra dimension. 

• They provide a realization to lowest order of the AdS/CFT correspondence. 

• They incorporate the self-gravity of the brane (via the brane tension) . 

• They lead to cosmological models whose background dynamics are completely understood 
and reproduce general relativity results with suitable restrictions on parameters. 

The review has highlighted both the successes and some remaining open problems of the RS 
models and their generalizations. The open problems stem from a common basic difficulty, i.e., 
understanding and solving for the gravitational interaction between the bulk and the brane (which 
is nonlocal from the brane viewpoint). The key open problems of relevance to astrophysics and 
cosmology are 

• to find the simplest realistic solution (or approximation to it) for an astrophysical black hole 
on the brane, and settle the questions about its staticity, Hawking radiation, and horizon; 
and 

• to develop realistic approximation schemes (building on recent work [4011 14311 13901 14021 14031 
[2471 [3Ml inn [201 [132]) and manageable numerical codes (building on [247l[364l [5n [202l [T37] 1 

to solve for the cosmological perturbations on all scales, to compute the CMB anisotropics 
and large-scale structure, and to impose observational constraints from high-precision data. 

The RS-typc models are the simplest brane-worlds with curved extra dimension that allow 
for a meaningful approach to astrophysics and cosmology. One also needs to consider gener- 
alizations that attempt to make these models more realistic, or that explore other aspects of 
higher-dimensional gravity which are not probed by these simple models. Two important types of 
generalization are the following: 

• The inclusion of dynamical interaction between the brane(s) and a bulk scalar field, so that 
the action is 



S^TTl [ d^xV^^J V'^^R ~ Kg^A^a-^* - 2A5($) 
2^5 J L 



branc(s) 



-A($) + + imattor 



(450) 
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(see [mi[2Tl l3Ml[T4l[m][T45] l4m [232l[TMl[TTn[^ 

13271 12221 13171 11501 [T5]). The scalar field could represent a bulk dilaton of the gravitational 
sector, or a modulus field encoding the dynamical influence on the effective 5D theory of an 
extra dimension other than the large fifth dimension [26] |98l EOT] [3641 [SD [55] [2431 (SHI |355l 

m\- 

For two-branc models, the brane separation introduces a new scalar degree of freedom, the 
radion. For general potentials of the scalar field which provide radion stabilization, 4D 
Einstein gravity is recovered at low energies on either brane (4111 13381 1285j . (By contrast, in 
the absence of a bulk scalar, low energy gravity is of Brans-Dicke type [156] .) In particular, 
such models will allow some fundamental problems to be addressed: 

— The hierarchy problem of particle physics. 

— An extra-dimensional mechanism for initiating inflation (or the hot radiation era with 
super-Hubble correlations) via brane interaction (building on the initial work in |134l 

\Tm{m[{m\m^\m[\m\\m\\m\ mm\m\m [Tm\m\ ). 

— An extra-dimensional explanation for the dark energy (and possibly also dark matter) 
puzzles: Could dark energy or late-time acceleration of the universe be a result of 
gravitational effects on the visible brane of the shadow brane, mediated by the bulk 
scalar field? 

• The addition of stringy and quantum corrections to the Einstein- Hilbert action, including 
the following: 

— Higher-order curvature invariants, which arise in the AdS/CFT correspondence as next- 
to-leading order corrections in the CFT. The Gauss-Bonnet combination in particular 
has unique properties in 5D, giving field equations which are second-order in the bulk 
metric (and linear in the second derivatives), and being ghost-free. The action is 



jABCD] 



2^5 J L V 

(451) 
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where a is the Gauss-Bonnet coupling constant, related to the string scale. The cos- 
mological dynamics of these brane- worlds is investigated in |122[ 13461 13481 13441 11621 1501 
[Ml 133 [3101 [M31 [m [m [m imi IZni [3II]. in HO] it is shown that the black string 
solution of the form of Equation (jl44p is ruled out by the Gauss-Bonnet term. In this 
sense, the Gauss-Bonnet correction removes an unstable and singular solution. 
In the early universe, the Gauss-Bonnet corrections to the Friedmann equation have 
the dominant form 

H'^ oc p^/s (452) 

at the highest energies. If the Gauss-Bonnet term is a small correction to the Einstein- 
Hilbert term, as may be expected if it is the first of a scries of higher-order corrections, 
then there will be a regime of RS-dominance as the energy drops, when cx . Finally 
at energies well below the brane tension, the general relativity behaviour is recovered. 

Quantum field theory corrections arising from the coupling between brane matter and 
bulk gravitons, leading to an induced 4D Ricci term in the brane action. The original 
induced gravity brane- world is the DGP model |132[ [95 ] [3471 1397] . which we investigated 
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in this review as an alternative to the RS-type models. Another viewpoint is to see the 
induced-gravity term in the action as a correction to the RS action: 



S 



2nl 



branc 



-g 
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(453) 



where /3 is a positive coupling constant. 

The cosmological models have been analyzed in [m] UMl [Ml USD El 1323 [3M1 ISZ3 
WM [71 ism UMl 13101 [Ml [MI]- (Brane- world black holes with induced gravity are 
investigated in |242J .) Unlike RS-type models, DGP models lead to 5D behaviour on 
large scales rather than small scales. Then on an FRW brane, the late-universe 5D 
behaviour of gravity can naturally produce a late-time acceleration, even without dark 
energy, although the self-accelerating models suffer from a ghost. Nevertheless, the 
DGP model is a critical example of modified gravity models in cosmology that act as 
alternatives to dark energy. 

The RS and DGP models are 5-dimensional phenomenological models, and so a key issue is 
how to realize such models in 10-dimensional string theory. Some progress has been made. 6- 
dimensional cascading brane-worlds are extensions of the DGP model. 10-dimensional type IIB 
supergravity solutions have been found with the warped geometry that generalizes the RS geometry. 
These models have also been important for building inflationary models in string theory, based on 
the motion of D3 branes in the warped throat [S31 1136) (see the reviews [2941 15^ and references 
therein). The action for D3 branes is described by the Dirac-Born-Infeld action and this gives the 
possibility of generating a large non-Gaussianity in the Cosmic Microwave Background temperature 
anisotropics, which can be tested in the future experiments [3981 1213] (see the reviews |5S1 1250] ). 

These models reply on the effective 4-dimensional approach to deal with extra dimensions. 
For example, the stabilization mechanism, which is necessary to fix moduli fields in string theory, 
exploits non-perturbative effects and they are often added in the 4-dimensional effective theory. 
Then it is not clear whether the resultant 4-dimcnsional effective theory is consistent with the 10- 
dimensional equations of motion \108i il09l i238i i251j . Recently there has been a new development 
and it has become possible to calculate all significant contributions to the D3 brane potential in 
the single coherent framework of 10-dimensional supergravity [31] [301 [Ml US] • This will provide us 
with a very interesting bridge between phenomenological brane-world models, where dynamics of 
higher-dimensional gravity is studied in detail, and string theory approaches, where 4D effective 
theory is intensively used. It is crucial to identify the higher-dimensional signature of the models 
in order to test a fundamental theory like string theory. 

In summary, brane-world gravity opens up exciting prospects for subjecting M theory ideas 
to the increasingly stringent tests provided by high-precision astronomical observations. At the 
same time, brane-world models provide a rich arena for probing the geometry and dynamics of the 
gravitational field and its interaction with matter. 
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